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Abstract. 

We present a complete derivation of the semiclassical limit of the coherent state 
propagator in one dimension, starting from path integrals in phase space. We 
show that the arbitrariness in the path integral representation, which follows from 
the overcompleteness of the coherent states, results in many different semiclassical 
limits. We explicitly derive two possible semiclassical formulae for the propagator, we 
suggest a third one, and we discuss their relationships. We also derive an initial 
value representation for the semiclassical propagator, based on an initial gaussian 
wavepacket. It turns out to be related to, but different from. Heller's thawed gaussian 
approximation. It is very different from the Herman-Kluk formula, which is not 
a correct semiclassical limit. We point out errors in two derivations of the latter. 
Finally we show how the semiclassical coherent state propagators lead to WKB-type 
quantization rules and to approximations for the Husimi distributions of stationary 
states. 
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1. Introduction 

Semiclassical approximations in phase space using coherent states have been discussed 
extensively for several decades. This attractive topic, a favorite of many theoretical 
physicists and chemists, turns out to be very difficult. In this contribution to its 
literature, we shall attempt to sort out and clarify the web of contradictions and 
inconsistencies that have characterized the recent state of the field. We shall do so 
for the simplest possible case, one-dimensional coordinate space, i. e. two-dimensional 
phase space. This is the case where it is relatively easy to check the semiclassical 
approximations. We have done work in higher dimensions as well, but we do not include 
it here, as it would only obscure the basic relationships and further lengthen the paper. 
The conclusions we have reached are stated in section 7, and the reader who is already 
familiar with the subject may jump to them now to get an overall view. Because the 
pitfalls are numerous, however, we shall follow a slower approach, a historical one in 
this introduction, and then a systematic and detailed one in the body of the paper. 

The study of semiclassical methods has two basic motivations. First, it provides 
approximations to quantum mechanical quantities in terms of classical ingredients. 
These approximations should be very good if the typical classical actions are much 
larger than Planck's constant. Interestingly, they are often fairly good even at very 
low quantum numbers. Second, semiclassical methods also help in understanding the 
quantum mechanical processes themselves, providing a more intuitive description. This 
description includes quantum mechanical interference, since both amplitudes and phases 
can be calculated semiclassically. 

The semiclassical approximation for the evolution operator, or propagator, in the 
coordinate representation has been known for more than 70 years and was first written 
by Van Vleck ||Van28|| . It is a complex number with a modulus and a phase. The main 



part of the phase is the action of a classical trajectory joining a given initial coordinate 
to a given final coordinate in a given time. Finding such trajectories is usually not 
a simple task. It is known as "the root search problem" and it gets more and more 
complicated as the number of dimensions increases. The modulus of the semiclassical 
propagator is related to the second derivative of the action with respect to these initial 
and final points. It measures the dispersal of nearby trajectories. Gutzwiller, among 
others, revisited this problem around 1970 [put71| , |Gut90|| , focusing on non-integrable 
systems and giving birth to the field of quantum chaos. Much progress has been made 



since then, particularly on the topological properties of Maslov indices ||Rob91|| and on 
the scars of periodic trajectories |[Hel84| , Pog88 |. 



But there is another interesting representation for the evolution operator, which 
seems at first sight to be more appropriate for comparisons with classical mechanics. 
This is the representation using the coherent states of a harmonic oscillator. They are 
gaussian states, localized in both coordinates and momenta, and therefore they can be 
thought of as quantum points in phase space. Although the exact coordinate propagator 
and the exact coherent state propagator are related by a simple change of representation, 
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the semiclassical approximations to them are quite different. One of the differences is 
that the classical hamiltonians with which the trajectories are calculated are different. 
Another is that the classical trajectories for the coherent state propagator are usually 
complex. Both semiclassical propagators involve trajectories with mixed initial and final 
conditions, hence both have the root search problem. 

This semiclassical coherent state propagator first appeared in the work of Klauder 
Kla78| , Kla79 , Kla87a | without a detailed derivation. Weissman Wei82 | extended the 



old semiclassical correspondence relations to the case of coherent state variables and 
presented a first derivation of the semiclassical propagator; his derivation was based 
on the general semiclassical machinery rather than on path integrals techniques. The 
possibility of a rigorous derivation using the latter is mentioned in several papers, but 
it does not seem to have actually been pubhshed, to our knowledge. The properties of 
the propagator were studied, however, for a number of fundamental quantum processes 



see, e.g., [pCav96a| , pCav96b| , pCav97| , |Gro98a|l ). A recent application to the semiclassical 



quantization of a system with classically mixed regular and chaotic dynamics [pch98| 
demonstrated the power of this approach. 

There is another important difference between the coordinate and coherent state 
representations. Due to the overcompleteness of the latter, the path integral for the 
coherent state propagator is not at all unique, and this non-unicity reflects itself in a 
large multiplicity of possible semiclassical propagators. For the quantum mechanical 
path integral, some good discussion of this variety was given in a review by Klauder and 
Skagerstam ||Kla85|| . Yet more ambiguousness arises because coordinates and momenta 
do not commute in quantum mechanics, and therefore there is no unique general way 
of associating a quantal hamiltonian with a classical one. The net result of all this 
is that many arbitrary decisions need to be made whenever one contemplates doing a 
semiclassical approximation. 

Without doubt, semiclassical approximations based on the propagator have been 
highly successful in chemical, molecular, atomic and nuclear physics. In spite of this, 
however, problems are turning up more and more often for which these methods are 
made very hard or inapplicable by the root search difficulty. These are usually problems 
in which the underlying classical dynamics is chaotic, which means that the number of 
contributing (real or complex) "root trajectories" can be extremely large. Consequently, 
people have attempted more and more to avoid mixed initial and final conditions. The 
ideal method is one in which one is given a coordinate and a momentum both at the 
initial time. Then the classical trajectory is unique, and the wave function evolves 
with time by following this unique thread. Such a method is called an "initial value 
representation" or IVR. Much work has been done on IVR's by many people, including 
Miller et al. ||Mil7CI| , |Sun97]] , Levit and Smilansky ||Lev77|| , Brumer et al. ||Cam92| , Pro95 



These IVR's do not involve coherent states. Klauder ||Kla87b|| also struggled with this 
problem with coherent states. Kay ||Kay94a] , |Kay94b| , [Kay97|| compared several IVR's 
numerically. A very popular IVR, based on coherent states, is the Herman-Kluk or HK 
propagator [|Her84 , |Her86| , Klu86 , Gro97 , Cjue98 , Gro98b |. For recent reviews see the 
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papers by Sepiilveda and Grossmann ||Sep96| , |Gro£ 



Finally we should note that the coherent state representation is not the only way 
to do quantum mechanics in phase space. A fascinating alternative is provided by the 
Wigner-Weyl representation |[Wig32 , Hil84 , Ber89|| . We shall not pursue this approach 
in this paper and we refer to a recent review article by Ozorio de Almeida |Pzo98|| . 

Contents of this paper : In section 2 we give a very complete derivation 
of the semiclassical propagator in the coherent state representation. The result is eq. 
( p.59| ) or eq. ( p.68|) . We have actually two different calculations of the path integral, a 
step-by-step calculation in section 2 and a more general method in the appendix. The 
latter is used again (in section 3) to perform a different integral. Section 3 discusses 
the variety of possible path integrals using coherent states and compares them. Section 

4 contains the derivation of our IVR, which is eq. ( |4.29| ). Our original purpose was to 
give a solid derivation of the HK formula but, when we were finished, we found a result 
very different from theirs, and in much better agreement with the expected behavior 
of such a formula. Also in section 4, we compare our IVR with Heller's IVR. Both 
have equal claims to being a correct semiclassical IVR, but they are different. Section 

5 returns to the HK propagator and points out the errors made in two papers where it 
was derived. It also explains why, in spite of being an incorrect semiclassical formula, 
HK still works (sometimes poorly) in some situations. In section 6 we Fourier transform 
the propagator from time to energy, which yields the Green's function in the coherent 
state representation. By looking for the poles of this Green's function, we obtain the 
quantization rule for the energy levels. By looking at the residues of the poles, we obtain 
approximate Husimi distributions for the stationary states. Finally section 7 contains a 
summary of our results and our conclusions. 
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2. The Semiclassical Coherent— State Propagator 

In classical mechanics, it is convenient to describe the time evolution by focusing on a 
trajectory in phase space. One candidate for a similar quantity in quantum mechanics 
is the operator which describes the time evolution in the coherent state representation, 
the coherent state propagator. In this section we shall construct the semiclassical limit 
of this propagator. The result is eq. ( |2.59| ). For convenience, we shall confine ourselves 
here to a single degree of freedom. The extension to higher dimensional systems will be 
discussed in the future. 

2.1. The Path Integral 

The coherent states \z) of a harmonic oscillator of mass m and frequency lo are defined 
by 

=e-^l"l'e"^^|0) (2.1) 
with |0) the harmonic oscillator ground state and 

In the above q, p, and are operators; q and p are real numbers; z is complex. The 
parameters 

h={h/muj)2 and c={hmuj)2 (2.3) 

define the length and momentum scales, respectively, and their product is h. We shall 
need the wave function of a coherent state in the position representation, which is 




exp ( —p{x — q/2) 



{x - qY 

262 J \h 

-}-{x/h-V2zf + '-{z-z* 



[2.A) 



Now we consider a system with Hamiltonian operator H{t). We restrict ourselves 
to "reasonable" Hamiltonians, i.e. we assume that H{t), written as a function of the 
creation and annihilator operators d^ and d, can be expanded into a power series of 
d^, d. The matrix elements of the evolution operator from time to time t in the basis 
of the coherent states ( |2.1| ) are 



K{z",t-z\Q) = (z"|Te-^/o^(*')di'|^/^ ^2.5) 
where T is the time-ordering operator. For time-independent H this is simply 

K{z\t-z',Q) = {z"\e-i'''\z'). (2.6) 
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In order to write the propagator as a path integral, we divide the time interval (0, t) 
into N small intervals of length t := t/N, and we write 



N-l 



K{z",t; z',0) = {zn 




H{t')dt' \ \zo) 



[2.7) 



where, for simplicity of notation, we identify \zj^) = \z") and |2;o) = \z'). If the time 
step r is small enough, the variable H{t') in the integral can be replaced by the constant 
H[tj), with some intermediate time tj G [jr, {j + 1)t]. Then the time-ordering operator 
becomes unnecessary and we get, with large A^, 



N-l 



j=0 



(2. 



Now we insert the unit operator, namely 
1 = 



IT 



.dxdy 
z) {z\ 



TT 



, , dz*dz . , 
\z) {z\ 



:2.9) 



everywhere between adjacent propagation steps. We denoted the real and imaginary 
parts of z by X and y, respectively. The sign of the last equation member on the right 
is actually undetermined, because the sign of the Jacobian depends on the order in 
which the variables are taken. To avoid any possible confusion, we state here that, in 
all integrations, dz*dz/2T[i actually means dxdy/TT. After the insertions, the propagator 
becomes a 2(A^ — l)-fold integral over the whole phase space 



K{z'\t;z\Q) 



N-l ,9 N-l 
d^Zj 



n 

N-l 

n 



TT 



ni( 

j=0 



,z,-+i|e ft 



dz*dzj 
2m 



(2.10) 



where / is defined by 



f{z\z) :=ln 



'N-l 
.3=0 



H(t,)T 



(2.11) 



with z := {zq, Zi, . . . , Zn) and its complex conjugate z* := {zq, z^, . . . , 2;^). The reason 
for writing z* separately from z as arguments of / is that they must be considered 
independent variables, because we are integrating over the two variables xj and yj for 
each j. Eventually, when we carry out the stationary exponent approximation, each of 
these two real variables will be allowed to become complex, which will result in four real 
variables for each j. Following Klauder [[Kla78|] , we transform the integrand e-^^^*'^-* as 
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follows 




{N~l 
j=0 



(2.12) 



with the abbreviation 



Using the coherent state overlap formula 



{z^+,\H{t,)\z,) 

{Zj+i\Zj) 



= exp + z*^^Zj - ^kjf I 



(2.13) 



(2.14) 



we write e 



as 



^ gxp ^ 

i=0 
Af-1 



1| |2 * -^1 |2 ^''"-17 



exp 



j=0 ^ 



(2.15) 



Later the limit N —>■ 00 (respectively r ^ 0) will be taken. Then the above 
summations will turn into integrals, and expressions (|2.8|) to ( p.l2|) would appear to 
be exact, were it not for the well-known problems attached to the meaning of such 
functional integrals. We see, however, that the treatment of time-dependent systems is 
(almost) identical to that of time-independent ones. 

2.2. The Stationary Exponent Approximation 



In the semiclassical limit of small h we can approximate the integral (|2.10|) , with eq. 
( p.l5| ), by looking for the places where the exponent / is stationary and replacing it in 
their vicinity by a quadratic form of its variables {z*,z). We call this the stationary 
exponent approximation or the gaussian approximation. In the literature it is often 
referred to as the stationary phase approximation or the steepest descent approximation. 
Strictly speaking, neither of these two names is quite correct. Our exponent / is complex, 
therefore it is not a phase. And "steepest descent" refers to a geometrical interpretation 
for the case of a single complex variable. 

Our approximation method involves going into the complex plane for the variables 
Xj and yj which are intrinsically real. Therefore it is important to keep clearly in 
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mind why we go through so many developments where we treat them as complex. The 
integrals we are after are over the real variables Xj , yj . Anything else we may do to 
calculate them is just mathematical tricks. The trick of going into the complex plane 
works only with analytic functions. Hence, when we have to do an integral over real Xj 
and yj, we shall first make sure that the function is analytic, and then we shall see what 
happens to this analytic function when we let the variables be complex. These words 
of caution may seem superfluous at this time, but they will turn out to be crucial later, 
when the need for further simplification arises. As we shall see in Sec. 3, confusion about 
this point seems to be what led Grossmann and Xavier |pro98b|| into error. As a start, 
one should test this analyticity requirement for the present integral ( 2.10| ). We already 
assumed in subsec. 2.1 that the operator H could be approximated by a polynomial in 
the operators and d. When we take its matrix element between the bra {zj+i\ and 
the ket \zj), we have to do an integral of the type 

{zj+i\H{tj)\zj) = J dx{zj+i\x) (Polynomial in x and ^) {x\zj) (2.16) 

where the two wave functions are given by (|2.4|) and its complex conjugate. The integral 
produces an analytic function of z*^-^ and Zj, with the additional factor 

1 



exp 



(2.17) 



which is the only place where the other two variables, Zj^i and z*, occur. According 
to (|2.13|) , this must then be divided by {zj^i\zj), which is given in (|2.14|) and which 



contains the same factor (|2.17|) , times another, never vanishing, analytic function of , 



and Zj. When the quotient is taken, the factor ( p.l7|) cancels out. Hence the "effective 
Hamiltonian function" 7ij+ij or T-C{zj_^_^, Zj;tj) of eq. (|2.13|) is an analytic function of 



the variables z*_^_i and zj separately, and it does not contain the other two variables Zj^i 
and z* at all ! 

Thus, the basic idea is to approximate the argument of the exponential in eq. 
( p.l5| ) by a second order Taylor expansion in the vicinity of the stationary trajectory. 
The resulting quadratic form in the exponent leads to a gaussian integral which can be 
done exactly. There may be more than one classical trajectory between the end points, 
each with its own contribution, which leads to a sum. For clarity during the derivation, 
this sum will not be written explicitly. We find the stationary points by requiring the 
vanishing of the derivatives of / with respect to z and z* separately, as mentioned earlier 

ir dHj+ij 



dl 

df 



dz. 



■i+i 



h dzj 
ir dHj+ij 
h dz*^^ 











3 



0, 



AT- 1 



Ar-2 . 



(2.18) 



We introduce new integration variables rj and rj* which describe the deviations from the 
points of stationary exponent: z — > z + r/ , — > + r^*, with the boundary conditions 



(2.19) 
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Now the exponent in eq. ( p. 15 ) 



j=0 

1 

2 



- -(4+1 + v]+i)izj+i + Vj+i - Zj - Vj) 



IT ^ 



(2.20) 



will be expanded into a Taylor series in {rj*,-!]) around the stationary points {z*,z) up 
to second order: 



1 



f{z* + 7]*,z + r])^ 2^ 1 2 (4+1 - 4)^i - 2^j+i{^J+i - zj) 



3=0 



IT 1 

^^j+i,i(4+i' ^i) + 2^4+1 - 4)^i 
+ ^(4+1 - 4)^i - ^4+1 - 



2-0 
1 



^j+i(^i+i - 



(2.21) 



•j+l 



^.+1 



+ 2(4+1 - 4)^i - 2^j+i(^^+^ - ^j) 



IT 

2h 



^3 + 



9z 



*2 ''i+l 

i+i 



f{z\z) + Y,{ 

j=0 



1 ffUj+i, 
^■+1 2 ^' n 



1 ir dHj+i 



- ::Vj+iV3+i + Vj+iV3 - T^VjV] 



4+1 9^j'4 



(2.22) 



2 

ir 

2h 



4 + 2 



•94+1% 



•94+1 



*2 '/i+l 



The terms of first order (second and third line of eq. ( p.22| )) vanish, when the boundary 
conditions ( |2.19| ) are taken into account, because of the stationary exponent conditions 
(|1|). Inserting eq. ( ^^221) into eq. (glOD yields, in view of eq. ( p.l9| ). 



' N-l 



dr]*dr]j I _ ^^^i+iJ ^2 ^ 9^^-1-3+1,3 

2h dz*l^ 



ir(z",t;y,0) = e/(^^^) / n 



i=i 



27ri 



exp 



2h dz^ 
3=0 ^ 



i,2 



^i+1 



VjV3 + 



V hdz*^^dzjj 



r^* r/, . (2.23) 
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2.3. Performing the Integrals 

At this point we start to carry out the integrations over the variables r]j,ri*. We have 
two ways of doing this. In the body of the paper, we perform successively // drjldrji, 
J f ^V2^V2, etc . . . , deriving eventually a recursion relation, which becomes a nonlinear 
differential equation when we go to the limit of continuous variables. In the appendix , 
we do it by writing the multiple integral in terms of the determinant of the quadratic 
form, we calculate this determinant with a pair of recursion relations, which turn into 
linear differential equations in the limit. The two methods are quite different, but they 
give the same result. The reader who wants to save time does not have to study both. 
In the appendix, then, we go on to use the same method to calculate a different path 
integral which arises in section 3. 

Here, we calculate the integrals J J drj*dr]j by applying the general formula 

'^^'^y ^Aix^+A2y'^+Aaxy+Bix+B2y ^ t q 4A^a^^ ^2 24) 



which is correct if the integrations are done along the two real they would be if 

X and y are the real and imaginary parts of r], respectively. Two comments need to be 
made here. First, the integral must converge. If we call — /ii and — /i2 the eigenvalues 
of the symmetric matrix of the quadratic form in the exponent on the lefthand side, 
convergence requires that /ii and /i2 both have a nonnegative real part. Second, the 
phase of the square root needs to be defined, since ^1,^2,^3 are complex numbers. 
This square root is also equal to the product ^/ziiy^, and the phase 0, of each ^/JTi 
must be chosen to satisfy —tt/A < (pi < tt/A. This phase rule is extremely important in 
determining the phase of the semiclassical propagator when one works in configuration 
space. For the present case of the semiclassical coherent-state propagator, it is less 
crucial because the phase can usually be determined by appealing to continuity in time. 

We must now rewrite this formula in terms of the variables rj and r]* rather than x 
and y. The transformation of variables is simple enough, but the paths of integration 
are totally changed, and the associated conditions need to be restated. To simplify 
notations a little, we call the variables u and v instead of 77 and 77*, respectively. Then 
the new formula is 



J„ J„. 1 +'>ia2+b2ai-bl''2'^3 

^^aiU''+a2V +a3UV + blU+b2V _ ^ a^-4aia2 ^2 25) 

27ri a/aI - 4aia2 

The integral is convergent if and only if the two numbers 



AH,2 = -03 ± 2^/0302 , 

which are the negatives of the eigenvalues of the symmetric matrix in the exponent on 
the lefthand side of eq. ( |2.24| ), not ( p.25| ), both have a nonnegative real part. The square 



root in eq. (|2.25|) is again the product ^/iiy^ and, once again, the phase (pi of each 



fj,i must be chosen to satisfy — 7r/4 < (pi < n/A. 
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For the rjl, rji integration, the parameters are 



2h dzl 2h dzf 

IT d^n2 



as 



1 



h dz^dzi 



V2 



and 



. 



(2.26) 



Convergence is assured since —as is real, positive, and much larger than oi and 02 (t is 
arbitrarily small). The correct branch of the square root is the one whose phase is close 
to 0. The result is 



K{zN,t;zo,0) 



l + 2i 



exp < 



h dzf 



'N-l 

n 

.i=2 



dr]*dr]i 
2'Ki 



h dztdzi 



X1V2 



*2 



1 + 2 



ir 8^2,1 
'h~dzf' 



Xi 



(2.27) 



j=2 3 



h dz*j^^dzj 



2h dzf 



The second set of integrals, / / d?72d?72, is done again with eq. (12^251) , but now a2 becomes 
more complicated and is expressed as a function of Xi. This kind of behavior continues 
at each stage of the integrations. If we call Xj the value of 02 when we do the integrals 
J J dr]jdr]j according to eq. (|2.25| ), then we find the following parameters 



ai 



i+i.i 



2h dz] 



as 



1 - 



h dzj^^dzj 



ij+i 



62 = (2.28) 



while the formula for 02 becomes a recursion relation for Xj 
Xq = (initial condition) 



h dzjdzj^i 



2h dzf 



1 



-Xj-\ 



J = l, 



AT- 1 



(2.29) 



Xj-i 



Once all integrations are done, the result is 

N~l 

K{zN,t;zo,0) = ef'^''''^l[ 



l + 2i- 



(2.30) 



h dz'j 



Xi 
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with Xj satisfying eq. (|2.29|) . Once again, the phase of each square root should be chosen 
close to 0. 

2.4- Continuous Variables 

The time has come to perform the limit N oo (respectively r — > 0). This gets rid of 
the approximations associated with the time discretization in eq. ( p.l5| ), and the discrete 
recursion formula (|2.29|) becomes a solvable differential equation. The stationary phase 
conditions ( p.l^ ) are in this limit identical to Hamilton's equations 

_idn 

h dz 
i dTL 

z = --— 2.31 



where 'H{z* , t) is the limit of T-^j+ij, eq. ( p.l3| ), and is simply given by 



n{z\z,t) = {z\H{t)\z) . (2.32) 

The question of boundary conditions presents us with a grave problem at this 
point. We have been approximating the propagator going in time t from coherent state 
{q',p') to coherent state {q",p"), and our approximation seems to involve the classical 
path between the two points in phase space. But there is no such classical path, real or 
complex, between these two points, generically speaking! The unique classical trajectory 
which goes through {q',p') at time does not in general go through {q",p") at time t. 
We have too many boundary conditions! 

One way out of this quandary was shown by ||Wei83|| . Actually, zn and Zq do not 



enter the equations of motion at all. Neither one occurs in eqs. (|2.15| ) or ( p.l8|) . Thus 



the problem is really to find a classical path going from zq to z*j^ in time t, and such a 
path does exist, but it is usually complex, in spite of the fact that (go)Po) and {qN,PN) 
are real. The path will be real only if (go,Po) and {qN,PN) happen to be on the same 
classical trajectory. Given that the intermediate values of q and p will now be both 
complex usually, it makes no sense to retain the notations z and z*, which can only lead 
to confusion. We replace these by the two complex variables u and v, which are now 
manifestly independent 

1 (q p 

Z ^ U = —= 7 + « - 

V2\h c 

.-.„^-L(|-,^), (.33, 



Note also the inverse formulae 



p=- — {u-v). (2.34) 
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The differential equations ( p.31| ) are now 

OH 

ihU = 

OV 

ihv = -— (2.35) 
ou 

with the boundary conditions 

m(0) = z' , v(t) = z"* , 

f (0) = nothing special ^ z'* , (2.36) 
u{t) = nothing special 7^ z" . 

v{0) and u{t) come out of the calculation and do not have any simple relation to z' 
and z", except in the special case when there exists a real trajectory going from z' to 
z" in time t. Then and only then do we get v{0) = z'* and u{t) = z" . Otherwise, 
the end points of the classical path are truly complex in phase space. To complete the 
definitions, we rename these end points as follows 

m(0) = u', t;(0) = v\ u(t) = u", v{t) = v". (2.37) 

Given the change of variables, we note the following differentiation rules, which follow 
from eqs. ( p.33|) and (|2.34| ) and will be needed later 



d dv d du d 1 f d d 

dq dq dv dq du -\/2 b\dv du 

d dv d du d i / d d . /'o qon 

dp dpdv dp du ^/2c\ dv duj 



— = —lb— + tc— 
dv a/2 V 9q dp 



d_ 

dv 
d 



, d . d 
b — — zc — 
du V dq dp 



(2.39) 



Let us rewrite f{z*, z), eq. ( p.l5| ), in the continuous limit and in terms of the new 
variables. Since / can also be written 



N-l 
j=0 



2t 



Z--Z* ^ 



Zj+i - z 



2r 



one might think that, in the limit t —>■ 0, this would reduce to 



dt' 



1 i 
-(?)m — iiv) — ■^'^(^) 



(2.40) 
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but that would be a mistake. Recalling the earher discussion, we realize that eq. (|2.4(]| ) 
assumes that u{t) = ujy, which is not equal to z", and it also assumes that f (0) = Vq, 
which is not equal to z'*. But in expression ( p.l5| ), un or zn was z", and vq or Zq was z'*. 
To correct this mistake, we must take out from the sum the two terms containing 
and fo, namely —^v^un and —^VqUq, and replace them by their correct values, namely 
— ^\z"f and — . Consequently, the value of / in the limit r should be 



/ 



dt' 



1 i 

-{vu — iiv) — -^'H{u, V, t') 



+ \{v"u" + v'u') - h\z"\^ + |zf ) . (2.41) 



Next we rewrite the product in eq. ( p.30| ), performing the limit N ^ oo and using 
ln(l +x) =x + 0{x^) : 



hm T] \ l + 2tl^xA = hm exp Vlnfl + 24^X, 



N-l 



(2.42) 



Altogether, the propagator is 



i^(z",t;/,0) = exp^ I dt'^(t')X(t') ^xp 



dt' 



- (vu — uv) — —H 
2 ^ ^ h 



+ ^iv'u' + v"u")-^{\zf + \z'f) 



(2.43) 



We still have to write the continuous form of the discrete recursion formula ( |2.29| ) 



for X{t). With 



1+bx 

'2h~dvf 



(2a + h)x + 0(x2) eq. (^7M) gives 



Xj-i 



+ 2 



In the limit N oo this leads to the nonlinear differential equation 

m- 

with the initial condition 



2hdv^ hdudv ^' hdu^ ^' 



(2.44) 



(2.45) 



X{Q) = . 



(2.46) 
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To solve this, we consider small variations in the solutions of Hamilton's equations eq. 
( p.35| ) around a given solution v{t), u{t): 



ou = — — — — — ou — — ^ „ ov 



Sv 



hdudv h dv"^ 
h du^ hdudv 

1 5u 



(2.47) 



It can be seen that the solution of eq. ( p.45| ) is X = • The time derivative is 

15u 1 5u 
2 5^ ~ 2 5^ 



X 



■ 6v 



2h dudv 6v 2h dv"^ 
i d^n i d^n 5u 



2h dv^ 



- 2- 



h dudv 5v 

I d^n 



I5u 

i d^n 

2h~du^ 



i_&^5u ^ i d'^n 
h du"^ 5v hdudv 

2 



5m Y 



(2.48) 



2h (9t>2 hdudv h du'^ 
which agrees with the differential equation ( p.45|) . The initial condition X(0) = can 
be satisfied by picking 5m(0) = 5u' = , 6v{0) = 6v' arbitrary. The integrand in the 
first exponential of eq. ( |2.43| ) can be transformed with the help of eq. (|2.47 ) 



i d^n 



X 



i d^nSu 



16v 
26v 



d^n 



h du"^ 2h du^ 5v 2 5v 2h dudv 
so that the first exponent of eq. ( |2.43D is 

,d^n 



-■^InSv 
2dt 



d^n 



2h dudv 



(2.49) 



exp 




du-' 



'it')X{t') 



exp 



5v' 
Sv" 




i d^n 

hdudv 



exp 



(2.50) 



It is understood that 6v'/6v" is calculated with the initial condition 6u' = 0. The phase 
of the square root evolves continuously with time, starting at for t = 0. Because all 
the numbers are generically complex, 5v"{t) does not usually have zeroes; therefore the 
phase is always well defined, barring an accident. In the end we obtain 



K{z",t;z',0) 



Sv' 



exp 



6v 



dt' 



-exp 




vu 



iiv) — —H 

h 



1 

2' 



(2.51) 



Just as a plausibility check, consider the limit t —>■ 0. In this case, the complex 

v" = z"*, and therefore K{z", 0; z', 0) = exp{0} exp{0 + 



u' 



Z' V 



path has u" 
\v'u' + \v\ 
overlap of two coherent states. 



llr'|2 
2 1^ I 



il r"|2 
2l^ I 



} = exp{- 



il r"|2 
2l^ I 



} = {z"\z'), which is the 
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2.5. The Complex Action 

In most discussions of hamiltonian mechanics, one gains much simphcity and under- 
standing by defining the 'action', which is the quantity entering in Hamilton's variational 
principle and in the Hamilton-Jacobi equation. In the present problem, this is true also. 
The action is complex in most cases, like the trajectories themselves and like the energy. 
It is given by the formula 



S{v" , u , t) : 



dt' 



ih 



{iiv — vu) — Tiiu, V, t') 



- j{u"v" + u'v') 



(2.52) 



The independent variables are v" and u', the two end variables which define the classical 
trajectory, and the time t ; u" and v' must be understood as functions of these three 
variables. Note that u and v are proportional to 1/ V^, since b and c in eq. ( p.33| ) are both 
proportional to ^/h. Therefore the only term in S which depends on h is 'H{u,v, t')dt'. 
We shall show in the following that this 5* is indeed the correct action for the boundary 
conditions we have. Given its definition, S should be an analytic function of its variables 
most of the time, since T-C{u,v) is analytic and the velocities it and v can be written 
as derivatives of Ti using Hamilton's equations ( |2.35| ). It would take an accident in 
the determination of the classical trajectory from the boundary conditions to produce 
a singularity. Similarly, the functions u" and v' should be analytic in v" and u' most of 
the time; in fact, according to eq. (|2.56|) , they are essentially the partial derivatives of 
S. 

Suppose that we make small variations 6v", 6u', 6t in each of the independent 
variables. This induces variations 6u", 6v' in u" and v'. It also induces variations 6u{t'), 
Sv{t') in the trajectory itself. The consequent variation in S is 



SS = J dt' 





— [VOU — UOV — VOU + UOV) 7—OV 7—OU 

2 ov ou 



ih 



+ 5t 



{v"5u" + u"5v" + v'5u' + u'5v') 
j{u"v" -v"u")-niu",v",t) 



We do two integrations by parts as follows 

/ dt'vSii = v"6u{t) — v'5u{0) — / dt!v8u 
Jo Jo 

- I dt'u6v = -u"6v{t) + u'6v{0) + [ dt'uSv 
Jo Jo 

Since the lower limit of the t'-integral is not changed, we have 

6u{0) = 6u\ 6v{0) = 6v'. 



(2.53) 



(2.54) 
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But the upper limit is changed from t to t + 6t, hence 

6u{t) = 6u" - u"6t, Sv{t) = 6v" - v"6t . 

Much simphfication occurs when we carry this back into eq. ( |2.53| ), and we find for the 
variation in S 

6S = dt' ' ^ \ e . ( ^ \ 







ihii — — ] Sv + { —ihv — — Su 
ov J \ ou J 

ih{u"6v" + v'6u') - n{u", v", t)5t . (2.55) 



The first fine says that paths satisfying Hamilton's equations ( p.35| ) have a stationary 
action S when the independent variables v",u',t are held fixed. The second line says 
that, for such a classical trajectory, the derivatives of S with respect to these variables 
are 

^ = -^^^"^ 1^ = -^^-'^ f = -nu", V", t) . (2.56) 

In the special case of a time-independent hamiltonian, T-C{u", v") is the constant energy 
£, and the last equality becomes 



We are now able to express most of the coherent-state propagator, eq. ( |2.51| ), in 
terms of the action function. By eq. ( |2.56| ) we have 



6v' i d^S 



5v" hdu'dv" 
Therefore the coherent-state propagator can be written 



(2.58) 



i^(.",M',o)=x:v^^exp|^l dt' (u^i 

exp|^^.(t;",«',t)-i(|;.'f + (2.59) 

The sum over u represents the sum over all (complex) classical trajectories satisfying 
the boundary conditions, since there may be more than one. We already mentioned (see 
before eq. ( p.l8| )) that we were going to suppress this sum for purposes of clarity, and 
we shall do so in the future again. 

There is one very important caveat. The hamiltonian function T-l{u,v,t) entering 
eqs. ( |2.52D and ( |2.59| ) is not the original classical hamiltonian of the problem, the quantal 
version of which we used when we first wrote the path integral ( |2.5| ). As eq. ( |2.32| ) shows, 
this "script" hamiltonian is a smoothed version of the original one, obtained by folding 
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it with a gaussian in phase space. This is one of the interesting features of formula 
( p.59| ). Another interesting feature is the first exponential, which we shall write 



efi 



where the quantity X, with dimensions of action, is defined by 

I{v",u',t) 



(2.60) 



(2.61) 



Though this X term was assuredly obtained in the past by various workers in the field, 
most authors seem to have been in a big hurry to forget its existence. One simple and 
compelling reason for not doing so is the following: when you calculate the harmonic 
oscillator (as we do in section 6), you get the exact answer if you keep X, but you don't 
if you drop it. 

There will be more discussion of X in sections 3 and 4. This, and only this, is 
the result of making the standard semiclassical approximations on the coherent state 
propagator ( p^.lO[ ). As far as we know, no formula equivalent to ( [^.59D has been in 
wide use before. Both Herman and Kay [Her86, Kay94a , Kay94b| ] claim, in words only, 
that a particular formula follows from a semiclassical treatment of the coherent state 
propagator. In both cases the claim is not justified and the formula is incorrect. In the 
next section, we shall see how other formulae can be derived in which the two special 
features above are either absent or reversed. In sections 4 and 5 we shall compare with 
two other formulae in the literature. 

For an additional note, we use eqs. ( p.34|) to write the integrand of X in terms of q 
and p 



dudv 



+ 



(2.62) 



2 dq^ 2 (9p2 

Given formulae (|2.3| ), we see that X is of order h. Whenever q and p are real, Ti is real, 
therefore d^H/dudv is real. Hence, if the classical trajectory happens to be real, X is 
real. 



2.6. The Tangent Matrix 

For the applications of our semiclassical formula ( |2.59| ) to be developed in sections 4 and 
5, we shall have to write the prefactor in terms of the elements of the tangent matrix 
M, which connects small displacements of the classical trajectory about the initial point 
at time zero to the evolved displacements at time t. Differentiating in eq. ( ^.56| ), but 
keeping the variable t constant, we can obtain the connection between initial and final 
displacements. In matrix form it is 



—ih 



/6u"\ 


(a a \ 




( 


[Sv'j 






[sv") 



(2.63) 
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with the notation 



du'dv" dv"dv" du'du' 

Solving for 6u" and 6v" in terms of 6u' and 6v' yields 





1 

1 




A 

\ 



A A \ 



-A. 



-It 



A 


— A 


\ 






/ 







(2.64) 



(2.65) 



If we call Muu,Muv, etc... the matrix elements of the tangent matrix (Note: unlike 
^uuy^uv, etc., these are not second derivatives) this last equation should also be 



(2.66) 



Therefore we have 



/ 6u"] 




(Muu 




l5u'\ 


ySv" ) 








\5v') 




i 




1 





h du'dv" M,„. 



and we rewrite the propagator (|2.59|) as 



^^e^vVj^S.{v\u\t)-'-{\zr + \zr) 



(2.67) 



(2.68) 



The square root in the prefactor has an undetermined sign. To know which sign is 
correct, one must start from t = when the square root is simply unity, and proceed by 
continuity along the trajectory, which can be done since M„„ never vanishes. We shall 
meet an important example of this in subsection 6.1. 

Another useful representation of the tangent matrix is in terms of the scaled 
variables q/h and p/c : 



5q"/h 
5p"/c 



^qq ^qp 



Tflpq TTLpp 



6q'/b 
6p' /c 



(2.69) 



The relation between the matrix elements above and those in eq. ( p.66| ) is: 



Muu = 


2 


{rUgg 


-\- Tfipp 


+ impg 


- imgp) 


Muu = 


1 
2 


[mgg 




+ impg 


+ irugp) 


Muu = 


1 








- irrigp) 


2 


[rrigg 




1771 pq 


Muu = 


1 
2 


{rUgg 




lj771jpq 


+ imqp) 



(2.70) 
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3. Ambiguities in the Choice of Path Integral 

There is more than one way of representing the propagator by a path integral in phase 
space in terms of coherent states. The way we chose in section 2 is just the one adopted 
by most workers in the field. Although each one of these different path integrals would 
give the same answer in an exact quantum mechanical calculation, they may differ 
when semiclassical approximations are made. In this section we shall discuss a number 
of such alternatives, and the implications of these ambiguities for the validity of the 
approximat ions . 

We begin with a qualitative remark. We introduced the path integral by splitting 
the propagator ( |2.5| ) into many segments, each with infinitesimal time-interval r, and 
then introducing the unit operator (|2.9|) between each pair of adjacent segments. But 
the basis \z) in terms of which the unit operator is written is vastly overcomplete (see 
Vou97|l for a recent discussion with references), and therefore there is an infinite number 
of ways of writing the unit operator in terms of the \z)^s. The way chosen for eq. (|2.9| ) 
is only one of them. We shall not pursue this approach to the ambiguous choice, but it 
does demonstrate the existence of an enormous arbitrariness. Instead, we shall discuss 
two other aspects of the problem. The first of these, through a different derivation of 
the path integral, leads to a different "effective hamiltonian" . The second deals with 
the arbitrariness in operator ordering when one goes from a classical hamiltonian to a 
quantum mechanical one. It will turn out that both of them will lead us to reconsider 
the significance of the first exponential (p.60|) in the semiclassical propagator ( |2.59|) . 



3.1. Alternative Forms of the Path Integral 



In their introduction to their overview of coherent states |[Kla85|| , Klauder and 



Skagerstam (KS) discuss two ways of arriving at a path integral, which they call "Path 
Integral - First Form" (p. 60) and "Path Integral - Second Form" (p. 69). The first 
form is the one we gave in subsection 2.1. The second form starts from the "diagonal 
representation" of the hamiltonian operator, namely 

H = J \z)h{z)^{z\ . (3.1) 

Given that we assumed (see second paragraph of subsec. 2.1) that H was either a 
polynomial in p and g or a converging sequence of such polynomials, this diagonal 
representation always exists. The notation h{z) for the function in eq. (|3.1|) is that 
of KS; we shall change it to H2{z*,z). This will be contrasted with the first-form 
hamiltonian function which we called 7i in eq. ( 2.32|) , and which in this section we are 



going to call Hi(z*,z), assuming no explicit time-dependence. Equation ( p.l|) is now 
rewritten 

H = j \z)H,{z\z)^{z\ . (3.2) 
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According to KS, the connection between Hi and H2 is 

Hi{z\z) = I f^\{z\z')fH2iz'\z') = J ^e-\^-^'\'H2iz'\z') (3.3) 

H2iz\ z) = (^exp Hiiz\ z) . (3.4) 

In addition to these two effective hamihonian functions, there is also the original classical 
hamiltonian Hc{z*,z). All three are different. 

A very common type of classical hamiltonian is Hc{q,p) = /2m + V{q) where, 
once again, we assume that the function V{q) can be approximated by a polynomial. For 
such a hamiltonian, the choice of quantum mechanical operator H is straightforward. 
Then we can work out the connection between He-, Hi, and H2 for simple monomials. If 
we call X and y the real and imaginary parts of z, respectively, this connection is given 
by the table 

H2 He Hi 



1 






1 


1 




X 






X 


X 




— 


1 
4 




X^ 


x^ + 


1 
4 


x^ — 


\x 




x^ 


x^ + 


3 ^ 
4*^ 


4 

X — 




3 

16 


4 

X 




3 


etc . 






etc . . . 


etc . 





(3.T) 
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with an identical table for monomials of y. We see that, up to cubic terms, we have 
He = ^{Hi + H2), but this does not remain true for higher powers. In a qualitative 
way, one can think of Hi as a smoothing of He, and of H2 as an unsmoothing. This 
idea becomes precise for monomials containing purely x or purely y. In the first case, 
we have 



Hcix) = y da;'e-2("-"')' 1/2(3;') 
Hi{x) = l/|y" dx'e-^^^-^'^"Hc{x') 



(3.5) 



showing that He is a smoothing of H2, and Hi a smoothing of He- The relations for 
pure functions of y are identical. If the monomial in He contains both x and y, then 
the problem of operator ordering arises; we shall discuss it in the next subsection. 

We shall now sketch the derivation of the semiclassical propagator for the second 
form of path integral. For simplicity, we exclude an explicit time-dependence in H. With 
fewer details, we repeat the steps of subsecs. 2.1, 2.2, and 2.3. We use the notation Ki 
for the propagator of the first form and K2 for that of the second form. The starting 
point is eq. ( p.8|) . Following KS (their eq. (6.11)), we write each of the N exponentials 
as 

-^^^ ^ f |2.)e-ir^^(-^^^)^^(z,-| . (3.6) 

vr 
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To facilitate the comparison between K2 and Ki, it is convenient to write e^'s^^* as the 
product of A^ — 1 factors hke ( |3.6|) rather than N. In other words, we take [N — l)r = t, 
which makes no difference in the hmit of large N and infinitesimal r. Then j in eq. 
(p.6|) goes from 1 to iV — 1. The complete propagator is 



K2{zn, t; zo, 0) = {zn\ {N — 1) factors similar to 

N-1 



3) ko) 



'j+i 



1%) • 



(3.7) 



j=0 



Contrast this with the first-form propagator, obtained from eqs. ( ^.101 ) and ( |2.12| ) 

K,iz^,t;zo,0) = / n ^ TT e-^^^(^^+-^^)(^,+i|^,) . (3.^ 



j=0 



Whereas in Ki the two arguments of Hi belong to two adjacent times in the mesh, 
the two arguments of H2 in K2 belong to the same time. This difference is important 
and results in different semiclassical propagators after one does the stationary exponent 
approximation. The stationarity conditions are found to be 

ir d 



h dzj 



H2{Z-,Z 



IT / * 

—Zj + Zj-i —Tr^H2[Z^ , Zj 



h dz* 











J 



3 



1,2, 



1,2, 



AT- 1 



AT- 1 



(3.9) 



Contrast this with the stationarity conditions for Ki, which are 

ir d 



— Zj + Zj_i — 



h dz. 
ir d 
h dz^- 



■Hiiz. 



Hi{z*,Zj.i] 











J 



1,2, 



J = 1,2, 



N-1 



N-1 



(3.10) 



The differences are subtle, but real. In the continuum limit, both sets of equations 
become the classical equations of motion, but for hamiltonian Hi in the case of Ki, and 
for hamiltonian H2 in the case of K2. In both cases, neither Zq nor zn appear in the 
equations, hence the trajectory is determined purely by zq and and is complex. The 
two actions Si and 5*2, coming from different hamiltonians, are different. 

Now we must calculate the prefactor. As in subsec. 2.2, the exponent in ( p.7| ) is 
taken at the points z* + 7]*, Zj + rjj and expanded to second order in the vicinity of the 
stationary points z*, Zj. The zeroth order gives an expression calculated on the classical 
trajectory, the first order vanishes by the stationarity conditions, and the second order 
is a quadratic form leading to a gaussian integration. The quadratic form for K2 is 



. 7V-1 



=1 >- 



1 2 



92 



1 .2 

h -ri 

dz*dz, 2'' 



dz*'^ 



N-l 



N-2 



(3.11) 
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According to eq. (|2.23|) the quadratic form for Ki is 

N-l 



ir X - 1 



^2 q2 



j^j- ^2 ^-2 

i=i j=i i=i 



(3.12) 



The comparison of K2 with i^i is easy if this gaussian integration is performed by the 
determinant method presented in the appendix of this paper. The details are given there 
and the result is the following. The semiclassical K2 is given by eq. ( |2.59| ) again, with 
two changes: (1) Of course, Hi is replaced by H2, and by S2; (2) The first exponent 
is — instead of +-|X . 

In those cases where the semiclassical approximation is expected to be good, one can 
then hope that the two changes cancel each other approximately. Both path integrals 
are exact originally and therefore they should give the same propagator. But one 
semiclassical propagator contains Hi and the other contains H2, which necessarily leads 
to different results unless the hamiltonians are purely quadratic. If both approximations 
are good, it must mean that the difference is cancelled by another one: the change in 
sign of the first exponent. 

Out of these considerations comes the justification for a procedure which has been 
used in some of the past literature: leave out the first exponential in ( p.59|) , or replace 



it by unity, and do the classical calculations using the original classical hamiltonian He 
instead of Hi (or 7i). This follows from the fact that, as eqs. ( |3.5| ) demonstrate. He 



can be thought of as more or less half-way between Hi and H2. It should then be 
associated with a first exponent which is also half-way between first-form and second- 
form, and that means zero. Another justification was given in [[Kur89|| , where it is 



shown that this is the way of getting the exact first-correction term in an expansion 
in powers of h. Hence the procedure is correct in the extreme semiclassical limit, the 
limit of high quantum numbers. It may not be as good for low energies, comparable 
with the energies of low excited states. We shall return to this question in a future 
publication. In our opinion, the best procedure for low energies is to take at face value 
the result ( p.59[) of the first-form path integral, using hamiltonian Hi and including 
the first exponential. This is because Hi is the smoothest of the three hamiltonians, 
and therefore the stationary exponent approximation has the best chance of being good 
in this case. There is one final question one might ask: besides the first form and the 
second form of path integrals, can one find a continuum of integrals which interpolate 
smoothly between these two? The answer is yes. Instead of the diagonal form (|3.6| ) for 
the infinitesimal propagator, one can write it in non-diagonal form thus 

-iHr f f I //\ / //I -^Ht\ /\ '^^^ / /I /o ION 

e " = \z ) {z |e \z ) {z \ . (3.13) 

././ vr 71 
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We can approximate the matrix element as in subsec. 2.1, which gives 

e-^^- ^ JJ |z")^^(/'|/)e-x^^i(-"*.-')^^ (^'1 . (3.14) 

After multiplying together many such expressions, one has to perform many integrals, 
twice as many as before. But half of the variables do not occur in the Hi functions, and 
the integrations over them are straightforward, the stationary exponent method being 
exact. This leads back to the first form of path integral. However, nothing prevents us 
from mixing the two recipes (|3.6|) and (p.l4|) in any proportions whatever, generating 
all possible interpolations between the first and the second form. We shall not pursue 
this approach further. 



3.2. Ambiguities in Operator Ordering 

In the above discussion, many of the statements were precise only if each monomial 
in the hamiltonian contained only q or only p. We shall now generalize to arbitrary 
monomials, which brings up the question of operator ordering. Given a classical 
Hc{q,p), there are many corresponding ff's, because q and p do not commute and 
their order matters. Hence there are many quantum mechanical problems and therefore 
many different results, the differences becoming smaller as h becomes smaller. The 
number of possible orderings, or combinations of orderings, is infinite, but three of them 
stand out. The normal ordering consists in rewriting Hq in terms of z and z*, replacing 
in each monomial z hj a and z* by (see eq. ( p.2|) ), and writing all creation operators 
to the left of all annihilation operators a, making each monomial look like CmnCf^'"^". 
The antinormal ordering does the same thing, but it writes the creation operators to the 
right of the annihilation operators, thus Cmnfl"^^'"- The third kind, called Weyl ordering 
or symmetric ordering, has several equivalent definitions. One says to write all possible 
orderings of the monomial, and then take the average of them all. Another definition of 
the Weyl operator Aw corresponding to a classical function A{q,p) is as follows. First, 
write A{q,p) double Fourier transform 

A{q,p) = jj dad/5 5(a,/3)e'(°'?+^P) . (3.15) 

Then the Weyl operator is obtained by changing q and p into operators in this formula 

Aw = II Mj3B{a,j3)S''^^^^^ . (3.16) 



Some easily-read references are |[Hil84|| , [[Kur89|| , and |[Vor89 



Thus we have three different ways of associating an operator with a function in 
classical phase space. We can call them An, A a, and Aw- Of the three. Aw is the most 
"reasonable" one in the classical limit. Conversely, there are three ways of associating 
a classical function with an arbitrary operator A, which are the inverses of the three 
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transformations above. The most interesting one (see [ ffil84 |) turns out to be the inverse 
of the Weyl transformation. It is 



s 



2 



Aciq,p) = / dse-^^Uq-- A q+-) . (3.17) 



This is called the Wigner transformation. Ac{q,p) is called the Weyl symbol of the 
operator A and is often denoted by Aw{q,p). The other two inverse transformations 
also have names. For the transformation which associates a normal-ordered operator 
to a classical function, the inverse transformation is called the Q-transformation, and 
the classical function is called the Q-symbol of the operator 

A(q,p) =^ normal ordering =^ Am 

. (3.18) 
AQ{q,p) Q-transformation <^= A . 

Similarly, the inverse of antinormal ordering is called the P-transformation, and the 
classical function is the P-symbol of the operator 

A(q, p) =^ antinormal ordering =^ A a 

. (3.19) 
Ap{q,p) P-transformation A . 

There is a very simple explicit expression for the Q-symbol 

AQ{q,p) = {z\A\z) . (3.20) 
For the P-symbol there is nothing as easy, but it is given implicitly by the requirement 

A = I \z)Ap{q,p)^{z\ . (3.21) 

Of the three symbols, the Q-symbol, which is analytic in both z and z*, is the smoothest. 
The P-symbol is the most likely to be singular. The Weyl symbol is in between. 

Coming back to our subject, we see now that Hi is the Q-symbol associated with 
the operator H, H2 is the P-symbol, and the function which we have called Hc{q,p) in 
the past should be the Weyl symbol Hw- As we had defined it for the pure monomials 
of q or p, it is indeed the Weyl symbol. And everything we have said about He so far 
will remain true in general, provided we define He as the Weyl symbol H^ of H. 

Let us summarize the results. We have three effective classical hamiltonians 
associated with the quantum mechanical H. Their smoothness decreases in the order 
Hi, Hw, H2. This is evident from the relations 



Hi{q,p) = l/f / dh'e-'\^-^'\'Hw{q\p') 
Hw{q,p) = Vf / d'^'e-'l^-^'l'i^2(g',pO • 



(3.22) 
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Going in the opposite direction, the relations are 



) 1 ( (3.23) 

In the extreme semiclassical hmit, i.e. for large quantum numbers, the best formula 
for the propagator is (|2.59|) calculated with Hw throughout and omitting the first 



exponential ( p.6(]| ). For low energy, the best formula is ( |2.59| ) as it stands, i.e. calculated 
with Hi and with the first exponential. The third formula, ( |2.59| ) calculated with H2 and 
with I replaced by —X (see after eq. ( p.l2|) and also the last sentence of the appendix). 



is valid too, though we do not know in what circumstances it might be expected to be 
better. All three propagators have claims to being called "the" semiclassical propagator 
in phase space, although there is no simple derivation of the propagator with the Weyl 
Hamiltonian from path integrals. They are all different since the classical hamiltonians 
are different. All three are exact for the harmonic oscillator and the free particle. But it 
would be wrong to calculate with Hi without including the I term, or with H2 without 
including —X, or with H\y including X. 

3.3. Powers of h 

In the following sections we shall use our semiclassical formula for the coherent state 
propagator in a number of different situations. In particular, we shall derive an Initial 
Value Representation for the propagator in section 4 and the Green's function in section 
6. In order to perform these calculations consistently, we must state precisely what is 
the 'philosophy' of our approximation, i.e., what terms must be kept and what terms can 
be discarded in calculations involving this semiclassical propagator. In order to do so, 
we have to understand the difference between calculations with the Weyl Hamiltonian, 
which does not involve h explicitly, and those with Hi = H or H2. The discussion boils 
down to understanding the semiclassical formulae in terms of powers of h. 

We start by recalling the nature of the stationary phase approximation (SPA) for 
a simple one-dimension integral. Let 

r+00 

A = / g{x) e^^(^) dx . (3.24) 



Assume that / has a single stationary point at X — Xq and that g{x) is a slowly varying 
function of x. The SPA can be applied if h is small. It amounts to expanding f{x) to 
second order about Xq while keeping g{x) constant in the vicinity of Xq. The integration 
is then reduced to that of a Gaussian, which is straightforward. The result is 



A-Ao:= J^^?(xo)e^+*^(-) (3.25) 
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where f^'^^ is the second derivative of / and s is its sign. This approximation neglects 
third and higher derivatives of /, as well as all derivatives of g. 

By keeping more derivatives, we can calculate the next non-zero contribution to 
the integral. This is done in appendix B. The result is that A can be written as 

A = Ao[l + ihR{xo) + 0{h'^)] (3.26) 



where -R(xo), given explicitly in eq. ( |B.10| ), involves the first and second derivatives of g 



and the third and fourth derivatives of /, all computed at If instead we were 

to write these corrections in the exponential, the result would look like this 



A 



V 1/^ ^1 

Therefore, since what we do all along is a quadratic exponent approximation, we do not 
have any hope of being able to calculate completely and accurately either the terms of 
order h in the integrals themselves, as ( ^.26 ) shows, or the terms of order /i^ in quantities 



occurring in an exponent with overall coefficient i/h, as shown in ( |3.27| ). However, since 
we actually do have terms of this sort in our results, more discussion needs to take place. 

The path integral calculation of the semiclassical propagator involves many 
integrals similar to ( ^.24[ ), which are evaluated by an appropriate quadratic exponent 



approximation. The role of f{x) is played by the action and that of Xq by the stationary 
trajectory. In the coordinate or momentum representations, neither the action nor the 
stationary trajectory depend on h. In these cases, as in the one- dimensional integral 
( p.24| ), corrections beyond the stationary phase approximation are of order h with respect 



to the semiclassical formula, as in eq. (|3.26|) . In the coherent state representation. 



however, both the action and the stationary trajectory depend on Planck's constant via 
Ti. Therefore, the quadratic approximation already involves 'nonclassical' /i-dependent 
terms, and an expansion of the propagator in powers of h similar to ( p. 261) becomes 
somewhat confusing. Interestingly, it is shown in appendix C that the non-classical 
terms in the action S cancel those in X up to first order in h. Therefore, the non-classical 
terms in the 'effective phase' S + 2 are of order h"^, which is beyond the precision of the 
approximation. It is not clear at this point whether these extra terms improve or not 
the semiclassical formula with respect to a pure 'Weyl' calculation. We shall come back 
to this point in section 6 when we discuss semiclassical quantization rules. 

In view of these results, this is what we should do in future calculations to be 
consistent with our semiclassical procedure. We must distinguish two parts to the 
integrand. One is the exponent that we make stationary. This is always multiplied by 
a factor i/h, and it usually consists of some kind of action. The other part, which is 
everything else, we call the prefactor. Then: 

(i) In the stationary exponent, we can drop terms containing third or higher order 
derivatives of ?i or S" with respect to u and v, since they have not been taken into 
account in the quadratic approximation. We should keep all terms of order h. They 
are the ones that are expected to cancel out. All terms of order can be discarded. 
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(ii) In the prefactor, we do not need to expand at all. Improving the result by expanding 
the prefactor is an illusion. In addition, the prefactor may have a phase, which of 
course we should keep. 



1 T7^Kv.n,„v.Tr onnc 



An Initial Value Representation 



30 



4. An Initial Value Representation 

The biggest difficulty in numerical applications using the coherent state propagator is 
the root-searching problem. This explains the recent popularity of the initial- value- 
representations (IVR), which avoid this problem. In this section we shall transform 
the propagator ( |2.59|) into a semiclassical IVR propagator with the same initial and 
final states as the Herman-Kluk (HK) propagator ||Her84|| . We had expected that the 



result would be the HK propagator itself, but it did not turn out that way: our IVR 
propagator differs radically from HK's. Along the way, we shall point out the mistake 
made by Grossmann and Xavier ||Gro98b|| in their attempt at a similar derivation. We 



shall return to the HK propagator in section 5, where we compare it to our IVR and 



point out that its original derivation in [|Her84|| also contains a mistake very similar to 
that made by Grossmann and Xavier. 

The basic purpose of an IVR formula is the following. We are given at time a 
wave function {x\ip{0)) in the usual configuration-space representation. We want to 
calculate its evolution, i.e. we want to know 

{x\m) = {^\Kit)\m) (4.1) 

again in configuration-space representation, K{t) being the usual Feynman propagator. 
We want to do this semiclassically, in terms of integrals over classical trajectories, and 
we want these trajectories to be specified purely in terms of the initial values of their 
coordinates. The alternative, mixed values of the coordinates, some being initial and 
some final, leads to unacceptable "root-search" difficulties. Among the many IVR 
formulae that have been proposed, it has been reported [|Kay94b|l that the Herman- 



Kluk expression ||Klu86|| is particularly easy to use and gives particularly good results. 



To use this formula, one must first transform the initial wave function to the coherent 
state representation (a.k.a. the Bargmann representation) by doing the integral 

(.1,(0)) ^ Ji.Wim) . (4.2) 

Then the propagation is carried out with a mixed propagator, which has coherent state 
coordinates initially and configuration coordinates finally, 

(xim) = J {x\K{t)\z')^{z'\m) ■ (4.3) 

The mixed propagator {x\K(t)\z') is the quantity calculated by HK. Our purpose is 
also to calculate it, but we obtain a different result (subsections 4.1 and 4.2). Our result 
is the one that follows naturally when one makes the type of semiclassical approximations 
that were made in section ^. Moreover, it has the essential property of conserving the 
normalization of the initial wave function, which the HK propagator does not have, as 
we shall see. Since the HK propagator has been so popular in the past, it is obviously 
desirable to make some numerical comparisons between the two formulae. In section 5 we 
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shall discuss the HK propagator in some detail and show some simple one-dimensional 
comparisons, for which the new formula is far superior. More tests are needed, including 
some for two-dimensional problems. 

It should be clear that the mixed propagator {x\K{t)\z') is simply a description of 
the time-evolution of a gaussian wave-packet. But a semiclassical approximation for 
this was proposed long ago by Heller ||Hel75|] . Obviously we need to compare it with 



our result. We do this in subsection 4.3. Once again the two formulae are different, but 
this time there is a strong resemblance, which we discuss in the light of the ambiguities 
encountered in section 3. Neither result is "better" than the other: they have different 
regions of validity. 

Finally, let us recall that the research on IVR formulae was originally motivated by 
the fact that the semiclassical propagator in coordinate space, long known as the Van 
Vleck formula, led to unpleasant root-search problems. In subsection 4.4 we examine 
the relationship between our and Heller's IVR with the Van Vleck propagator. 

4.I. A mixed representation 



Our starting point is the coherent state propagator of eq. ( |2.59| ). Our expression for the 
mixed propagator is then 

{x\K{t)\z') = J{x\z")'^{z"\K{t)\z') 

^ 7r-4r2exp<^-^— + x-^U (4.4) 



e 



27rn M 262 n \ 2 



i d^S 



h du'dv" 



u'=z\v"=z"- 



exp|^S(t;",.',t)-i(|.'f + 



Before one does the integration over dp"dg", it is good to remind oneself of the philosophy 
spelled out in the second paragraph of subsection p.2| . The integrals to be done are 
actually two real integrals over dp" and dg", and if the functions are analytic in p" and 
g", then it is all right to deform the contour into the 4-dimensional space of complex p" 
and complex q" . But there are still two integrals and, therefore, in this 4-dimensional 
space the integration runs over a 2-dimensional surface. Applying these thoughts to the 
last member of eq. ([4.4|) , we see that the argument of the exponential in the first line 
is obviously analytic, and we can rewrite it in terms of u" and v" without problem. In 
the second line we have 



2 V 62 c 



(4.5) 



which is analytic and can also be written u"v" . Also in the second line, we have the 
function S{v" , z' , t) and its second derivative. For them, the argument of analyticity 
and change of variables was made at an earlier stage already, when the coherent state 
propagator was derived in section 0. Recall that S{v",u',t) is a complex function 
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of the real variables p",q",p',q' which depends on these variables only through the 
combinations v" and u'. The fact that it is an action calculated along a complex 
trajectory is really not relevant to doing the integrals. The thing which is relevant 
is that it is analytic in v" and u', and therefore, if we have to do integrals over p", q", 
p', or q', we may continue them into the complex planes of these variables. Everything 
that was just said about S{v",u',t) can be repeated word for word about I{v" ,u' ,t), 
which also occurs in the last line of (|4.4|) . Finally, another very relevant thing is that 
we may wish to talk about the derivatives of S, and therefore we give them names as 
follows: 

Then, in terms of m" and f which are more convenient variables than p" and g", here 
are the integrals we have to do 



{x\K{t)\z ) = / — TT ib 2 Qh-^^^ - 



2ni \ h du'dv" 



ui=z>v"=z"* 



X 



2 



exp^ + ^xu" - \u"' - u"v" + '-S{v\ z\ t) - i|zf 



(4.7) 



We are going to do these integrals by the stationary exponent approximation. First, 
we must look for the stationary point of the exponent. We do not include the X term in 
the exponent, because that would involve calculating third order derivatives of 7i (see 
section 3.3). To find a point in 4 dimensions, one needs 4 real equations, or 2 complex 
ones. If we call F the exponent in the second line of eq. ( [4.7|) , these two complex 
equations are 

— = -^x - u" -v" = . 
ou" b 

There are two crucial comments to be made here. One is that S does not depend on u" 
at all, and therefore the second equation does not contain any derivative of 5*. In fact 
the second equation, according to eq. ( |2.33| ), says very simply 

q" = x . (4.9) 

The second crucial comment is that u" and U" are not the same. One is the independent 
variable u" . The other, U", is a function of the other independent variable v" . The 
fact that U" is obtained by calculating a certain complex trajectory is interesting, but 
irrelevant to the integration problem: it is some function of v" . On the other hand, 
u" is an independent variable; it can be anything, irrespective of the values of v" and 
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u' . Therefore the first equation, taken by itself, defines a 2-dimensional surface in 
4-dimensionaI complex space. One of the points on this surface, and only one, is the 
end point of the real trajectory which begins at p', q' . To find the stationary point, one 
must combine this equation with the second one. These two important facts are the 



ones that were missed by Grossmann and Xavier ||Gro98b|| . They totally ignored the 
existence of the second equation. And they claimed that the first equation said that 
the end point of the real trajectory was the stationary point. This is not so, unless one 
happens to pick x equal to the value of q for this end point. For all other choices of x, 
the stationary point is the end point of a complex trajectory. Its q" is real and equal to 



X . But its p" is complex. We give its value in eq. (|4.13|) 



We could take this result and use the stationary exponent approximation in the 
vicinity of the complex stationary trajectory. There are two defects to this approach. 
One is that finding this complex trajectory is once again a root-search problem, since 
it is specified by one initial coordinate z' and one final coordinate x . The other defect 
is the inconvenience of having to do classical mechanics in the complex domain. It 
is obviously preferable, if one can, to set up the practical applications in a way that 
involves the actual calculation only of real trajectories. This can be done, and at the 
same time the wished-for IVR character is recovered. The idea is that the complex 
trajectories which are in the vicinity of a real trajectory are a little bit like the gaussian 
which is in the vicinity of a stationary point when you do stationary phase or steepest 
descent integration. In the latter case, most of the contribution to the integral comes 
from the vicinity of the stationary point, so that you can approximate the exponent in 
this vicinity by a quadratic, even though it is not really a quadratic. You do a Taylor 
expansion of the exponent near the stationary point and you keep only the first two terms 
(actually, the first term vanishes). Similarly, when summing over the whole bunch of 
classical trajectories, we expect most of the contribution to the integral to come from 
the vicinity of the real trajectory, so that we can approximate the complex trajectories 
nearby by doing some sort of Taylor expansion to second order. The assumption is that 
the contribution of a complex trajectory falls off gaussian-like as it gets farther away 
from the real trajectory. This assumption is correct at least for simple systems like the 
free particle and the harmonic oscillator and it was shown numerically to hold for the 
quartic oscillator as well [|Xav96a . 



Returning to eq. (|4.71 ), we see that three terms in it prevent us from doing the 
integral exactly: the two terms containing the function S, one in the exponential and 
one in the prefactor, and the term containing X. In accordance with the ideas expressed 
above, the 5* in the exponential, which is the value of the action for the complex classical 
trajectory, will be expanded to second order in powers of v" — Vr, where Vr (or Zr*) refers 
to the final point of the real classical trajectory issued from z', thus 

^5'(f", z',t) ^ ■^'^('^'■' -2'' + "^riv" - Vr) + ^li'v" - 'Vrf (4-10) 
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with 



Ur 



i OS 
hdv^ 



(4.11) 



and 



7=1: 



h dv 



Jl2 



qp 



''pq 



m. 



■pp 



rugq + inia^ - im^a + m. 



(4.12) 



''qp 



pq 



''PP 



Here we have used eqs. ( |2.65[ ) and ( p.70| ) relating second derivatives of the action and 
elements of the tangent matrix. A similar expansion to first order in v" — Vr can be done 
in the first equation ( [4.8| ), which allows us to calculate p", the momentum at the end 
point of the complex trajectory 



P 



Pr + i 



7 c , 

- [X 



Qr 



(4.13) 



1 + 7 & 

As for the second derivative of S under the square root in the prefactor, it presents 
us with a bit of a problem. We know it (with its factor i/h) to be equal to the inverse 
of My^, taken at u'{= z') and v". If we were to expand this in the vicinity of the real 
trajectory, we would be taking derivatives of the tangent matrix, which is itself made 
up of second derivatives of S, and going to higher order than anyone ever goes in this 
kind of semiclassical argument. Common practice would say: just replace it by its value 
at the stationary point. Unfortunately, as we already know, the stationary point is not 
the real trajectory. It is a complex trajectory, hopefully rather close to the real one but 
not simple. Moreover, if we tried to do that, we would be back into the pitfall of mixed 
initial and final conditions. Hence the only reasonable thing to do for this presumably 
weakly-varying term is to use its value for the real trajectory, which is the place that will 
turn out to give the maximum contribution anyway. (Of course there are phase-space 
points q" and p" for which the difference v" — Vr is not small. But at these points the 
propagator itself should be negligible). The very same argument holds for the X term, 
since it already involves second derivatives of Ti. Hence we shall replace it by the value 
it takes for the real classical trajectory issued from {q',p'), and we shall call it X,.. 

Here is the complete formula now, in terms of the variables u" and w" := v" — 



{x\m\z') 

(iu"(kv" 



TT ih 2eri ''M^^ 2 exp 



2m 



1 //2 1 //2 

■ exp "2"" ~'~ 2^"^ 



u"w" + 



V2 



-X 



Vr 



U" + UrW" 



(4.14) 



where it is understood that M^^ is taken for the real trajectory. The integral is now a 
pure gaussian and we can do it using formula (|2.25|) , with 



ai = -- 



a2 



;7 



as 



-1 



Aaia2 



1 + 7 • 



V2 



-X — Vr 



Ur 



(4.15) 
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The two quantities 

Ati,2 = -03 ± 2y/aia2 = 1 ± 2^/7 
must have a non-negative real part, which requires 

\lm^\ < 1 . 



(4.16) 



(4.17) 



We have shown, with a fair amount of algebra, that this condition is indeed satisfied. 
According to eq. ( |2.25|) , the value of the gaussian integral, i.e. the second line of eq. 
( TO , is 



1 



1 



exp 



vr+7 1+7 



2 +^r[—X-Vr +-\—X~Vr 



(4.18) 



All we need to do now is to simplify the result. 

The new prefactor (1 + 7)"^^^^, multiplied together with the other prefactor 
{Myy)~^/'^, yields {M^y + M„t,)~2, which is the same as (m^g + imqp)~^ . Now let us 
rewrite S{z*,z',t) in terms of the usual action of Hamilton. Equation ( p.52|) says 



dt' 



1 i 

-(?)u — iiv) — ^5 



(4.19) 



We rewrite the first term in the bracket in terms of p and q: 

complex conjugate 



[vu — uv) 



qp qp 



{qp - pq) 



he be/ h 
since be = h . With an integration by parts we have 



* ft t 

dt'{qp — pq) = 2 dt' qp — pq 
Jo 

= 2 pdq - prqr + p'q' 
Jo 



(4.20) 



(4.21) 



where and Pr are the final points of the real trajectory starting at {q',p'). As for the 
last term ^l-ZrP of eq. ( f4.19| ), it will be convenient to write it ^UrVr ■ Expression ( |4.19| ) 
is therefore 



-S{z,,z,t) - -\z I 



« / , 1 



'^Sh - ^^Prqr + ^^pq + ^U.Vr 



where Sh is Hamilton's action 



Sh= !\pdq-ndt') 

J i 



(4.22) 



(4.23) 
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for the real trajectory. Carrying this into eq. ( 4.14 ), whose second hne is eq. ( [4.18 ), we 
find for the complete exponent in {x\K{t)\z'), in addition to the term: 



1 



1 + 7 



, 7 2 , V2 / ^ , 

- UrVr + -V^ + —x{Ur - 'JVr) + 'J-^ 



(4.24) 



This is where the serious work of simplification begins. 

We start by gathering all the terms containing x and we complete the square in x, 
which gives 



1 — 7 X 2{Ur — ■yVr) 



+ 



262 1 + 7 bV2 1+7 



1 — 7 / X Ur — 7f r 



+ 



(4.25) 



1 + 7 \hV2 1 - 7 / 1-7^ 
The first term on the right hand side of eq. ( [4.25| ), written in terms of and g^, becomes 



1 — 7 / X 



qr . 1 + 7 Pr- 

— I 



1 + 7 \bV2 bV2 1-7 cV2 

11-7/ x-qrV i 1 1 + 7 



(4.26) 



Return now to expression ( [4.24| ) and gather all terms quadratic in {ur,Vr), or in {qr,Pr), 
including the last term on the right hand side of ( 4.25| ): 



1+7 ,72, iUr~-fVry 
UrVr H UrVr H H 

2 2 2'^ 1-7 



(4.27) 



When this is written solely in terms of and Pr, much simplification occurs and one is 
left with the single term 



1 1 + 7 i4 

'2 1 - 7 C2 



(4.28) 



which cancels the last term of (|4.26| ). 

This is the end of the simplifications. The final formula is 



{x\K{t)\z') 



TT 40 2 



im. 



e'i'^'' exp 



qp 



1 ^ — 1 f X — Qr 



2 1+7 V b 



+ ^lpr{x-qr) + ^p'q' + Sh 



(4.29) 



This formula, and not the HK formula, is the logical consequence of transforming 
the semiclassical coherent state propagator to a mixed representation by applying the 
standard semiclassical approximations. 
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4-2. Some properties of the mixed propagator 



According to eq. ([4.29 ) the mean position of the packet (g) and its mean momentum 



{p) are given by qr and pr, the real classical trajectory originating at q',p'. Its squared 
width {{q — qrY) comes out of the real part of the coefficient of the gaussian, which 
should be and is negative. Since we have 



1-7 _ f^pv - if^pq _ 1 - i{mppmqp + mqgmpg) 
1 + 7 rnqq + iuiqp m^g + m^^ 



(4.30) 



we find 



52 

{{q - Qrf) = y [nilq + m^J = Agl,, (4.31) 

where Aq'^^^^^ is the classical spreading of a gaussian initial ensemble corresponding to 
the initial phase space distribution \{p,q\p',q')\'^- In similar fashion we have 

HP - Pr?) = J Hq + mlp) = ApL,, . (4.32) 



We see incidentally that the denominator of the prefactor in ( |4.29| ) can never vanish, 
and therefore the mixed propagator can never be singular, because the determinant of 
the tangent matrix must equal unity. This follows from symplecticity, which in one 
dimension is the same as Liouville's theorem. 

Another straightforward calculation yields the normalization of the packet 



N{t) = / dx\{x\K{t)\z')\ 



2 



Im? + ml [dxexp<-l 



V <iPj—-^]^2 



1-7 1-7* 
+ 



1 + 7 1 + 7^ 



{x - qrf 

62 



(4.33) 



Normalization is conserved and equal to unity at all times, which is as it should be. 

It may not be superfluous to mention once again that the sign of the square root 
in eq. ( [4.29|) is to be determined by continuous displacement along the trajectory, given 



that this square root is unity at t = 0. 



4-3. Comparison with Heller's IVR 



Now we proceed to the comparison with Heller's approximation mentioned earlier. 
This approximation is also known as "the thawed gaussian approximation" or TGA. 
A different derivation of the TGA was given later by Kay [ [Kay94a| | and used by him in 
numerical comparisons with other approximations [ [Kay94b| |. Although Heller does not 



give a final formula, while Kay does, we shall stick to Heller's presentation, as we do 
not find Kay's arguments convincing: they seem to be based on convenience and (very 
limited) numerical agreement, rather than solid basic principles. Heller's paper assumes 
a hamiltonian of the form p^ /2m + V{q). Here we shall consider a general hamiltonian 
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H{p, q) and we shall carry the derivation all the way to an explicit formula, which Heller 
does not do. The result will be that Heller's approximation leads to a formula identical 
to our eq. (^4.29|) except for two differences: (1) the classical hamiltonian that must be 
used to compute the trajectories, instead of being the smoothed H or the Q-symbol 
defined in section 3, is the Weyl symbol Hw of the quantum mechanical operator; (2) 
the term e^'^'' is absent. Hence the discussion of section 3 returns to the fore: there exist 
indeed different approximations, all legitimately derived, and the question becomes one 
of deciding under what circumstances one or the other can be expected to be better. 
We already saw in section 3 that, in the extreme semiclassical limit, when h becomes 
very small. Heller's approximation is expected to be best | Kur89|| . At low energy, on 



the other hand, we expect (|4.29| ) to be best. We gave a qualitative argument for this in 
section 3, but we have also performed several numerical comparisons, which we reserve 
for a future publication, as this paper is probably too long already. 

Heller's idea is to assume that the original wave-packet {x\z') remains gaussian as 
it propagates, but that the parameters of the gaussian change with time thus 

{x\K{t)\z') = exp [ait) (x - q{t)f + p{t) {x - q{t)) + (3{t)] ■ (4.34) 

He further assumes that q(t) and pit), which are the expectation values of the operators 
q and p, follow a real classical trajectory for some classical hamiltonian H{q,p). On the 
other hand he takes a{t) and to be complex, so that there are 6 undetermined 
real functions in the formula. The problem is to determine them. The exact 
quantum mechanical packet obeys the Schrodinger equation with a certain quantum 
hamiltonian operator H{q,p). What should be the connection between H and H7 The 
most straightforward assumption, and the one that works for the simple problems of 
elementary quantum mechanics, is to say that they are Weyl-related, i.e. H is the Weyl 
operator associated with H, and H is the Weyl symbol associated with H. But we do 
not want to make this assumption explicitly: it should come out of the calculation, and 
it will. 

Heller argues that, since the wave-packet is small, the only part of the quantum 
hamiltonian that matters is the part which, in phase space, refers to the vicinity of 
the classical trajectory, i.e. the region where the wave-packet is appreciably different 
from 0. In that small region it is permissible to expand the quantal hamiltonian up 
to second order in the variables q — q and p — p. This "local quadratic expansion" is 
the embodiment in this case of the semiclassical approximation, quite similar to the 
stationary exponent approximation in section 2 and to the second order expansion in 
eq. ( ^4.101 ). Hence we expand the hamiltonian operator as follows 

H{q,p) = H{q,p) 

+ H,{q-q) + H,{p-p) (4.35) 
+ i {H^^{q - q)^ + Hpp{p - p)^ + H^^ [{q - q){p - p) + {p - p){q - q)]} 
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where q is the operator 'multiphcation by x\ p is —ihd/dx, and Hg, Hp, Hgg, Hpp, Hgp 
are the first and second derivatives of H{q,p), which are functions of q and p and do not 
contain any operators. The crucial point to note here is that, if formula (|4.35|) is true, 



the Weyl symbol of the operator H{q,p) is simply H{q,p), without any additional terms. 
This comes out from repeated applications of the Wigner transformation (|3.17|) , where 
one does the integrals by using the standard properties of the Dirac delta-function and 
its derivatives. Denoting the Weyl symbol with the subscript W, one finds 

{q)w = q {p)w = p 

ie)w = q' {f)w=p' (4.36) 

{qp)w = qp+ ^ih {pq)w = pq - ^HT' 

etc .... 

Note that the Q-symbol of the right hand side of ( [4.35| ) would be very different. For 
instance, while the Weyl symbol of (g — q)^ vanishes, its Q-symbol is the squared width 
of the packet. Thus we have chosen H{q,p) to be the Weyl symbol of H, and this 
is not a trivial choice. It is however a purely arbitrary choice for now. The crucial 
moment will come when we prove that q{t) and p{t) obey Hamilton's equations for this 
particular classical H{q,p). This moment is close at hand. Meanwhile, we note that the 
last bracket [■ ■ ■ ] of ( [4.35| ) can be written 

— q){p — p) — ih . 

Now we try to satisfy the Schrodinger equation ihdtp/dt = Hip. For ?/'(t) we 
substitute the wave-packet ( [4.34|) and for H we substitute ( [4.35| ). We find 



1 dip 

—ih—- = —a(x — qY + 2aq{x — q) — pix — q) + pq — $ 
ip ot 

^HiP = H{q,p) + Hg{x -q) + ^Hgg{x - qf + 2Hpa{x - q) (4.37) 



+ Hpp \2a^{x - qf - iha] + Hqp 



2a{x — qf — -ih 



Matching powers of {x — q), we get the following 3 complex equations to determine the 
6 unknown real functions 

a = -2a^Hpp - 2aHgp - ^Hgg (4.38) 
2aq-p = 2aHp + Hg (4.39) 
$ = pq- H{q, p) + ihaHpp + ]^ihHgp . (4.40) 

We start with the imaginary part of (|4.39|) , which gives 

q = Hp. (4.41) 
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Then the real part of the same equation says 



P 



(4.42) 



Hence q{t) and p{t) follow a classical trajectory of H ^ provided that the latter is the 
Weyl symbol of H. Two equations remain, ( [4.38 ) and ( 4.40|) . 

To solve ( |4.38| ) we notice that this equation becomes identical to that for eq. 
( p.45| ), if we identify v with g, u with p and jHi, with H. Following the calculation that 
leads to eq. ( p.48|) , we find immediately 



a 



U_p 
26q 



where Sq and Sp are deviations from the classical trajectory, which satisfy 



(4.43) 



5q = Hqp5q + Hpp5p 
5p = -Hqg5q - HpqSp 



Using the tangent matrix eq. ( p.69| ) we get 

_ 1 i'rnpg5q' + rUppSp' 
2 mgg6q' + \mqp5p' 



1 Irupq + 2mppa' 

2 rugq + 2\mqpa' 



where 



a 



1 ^r)' 

J- = th/2b^ = ic/2h . 

^ (J 



Therefore 



a = — — — — 

2b + iniqp 



2b niqq + iiriqp 



ic 1 — 7 



2b 1 + 7 ■ 

Equation (^4.40|) can now be integrated. It can be written 

I3=pq- H{q,p) + YJ^^PP + y^9P 

where L is the Lagrangian. Integrating both sides gives 

P = Sn + jhgSq 



(4.44) 
(4.45) 



(4.46) 



(4.47) 



(4.48) 



(4.49) 



(4.50) 
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up to an additive constant, Sh being Hamilton's action (for the Weyl hamiltonian) , and 
with 



6q = niqqSq + \nqp5p' = 5q'{mqq + —mqpa') . 
The initial value 5q' is determined by /3(0) = q'p' /'i = ih/2logSq' , which gives 



6q' 



-iq'p' /h 



Substituting eqs. (|]Mp, ( ^3^ , and ( ^TTj ) into (|CT| ) gives 



(3 = SH + po^ 
= Sh + Y 



'^'H + ^ + y log (mgg + «mqp) 



or 



a. 
e 1 



A/mg,^ + im. 



,j-{SH+p'q'/2) 



qp 



(4.51) 



(4.52) 



(4.53) 



(4.54) 



up to a multiplicative constant. 

With these expressions for a and /3, one sees immediately that the Heller wave- 
packet ( 4.34 ) is identical to ours ( 4.29| ), except for the different hamiltonian (Weyl's) 



and the absence of the exponential. Equation ( |4.33| ) shows that the normalization of 
the Heller packet is conserved at all times. 



4-4- Recovering Van Vleck's Formula from the IVR 

The semiclassical limit of the evolution operator in the coordinate representation is 
given by the well known Van Vleck formula |[Van28|| . It is desirable, therefore, that other 
semiclassical representations of this operator reduce to Van Vleck's when transformed 
back to coordinates. In this subsection we shall calculate 



K{x",t;x',0) := {x"\K{t)\x') = I {x"\K{t)\z')^{z'\x') 

71 



(4.55) 



for both mixed propagators. Heller's and ours. The integration over q' and p' will 
be performed by the stationary exponent approximation. We shall see that Heller's 
approximation recovers Van Vleck's formula exactly. Our semiclassical approximation 
recovers it only in the limit of small h (see discussion in section 3), but this is also the 
only limit in which the Van Vleck approximation can be justified. 
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We start by inserting the Heller propagator in eq. ( 4.55| ). Using expression ( p.4|) for 

{x\z) we get 



-^Hciicr(a;",)f:;x',0) 



dq'dp' vr 26 1 



X exp 



X exp 



262 



p{x' — q /2) 



(4.56) 



We call ^ = ^{q',p') the exponent in the second and third lines of ( [4.56|) . The stationary 
conditions are d^/dq' = and d^/dp' = 0. We write ^ explicitly first: 



^ 262 



1-7 
1 + 7 



/\2 



{x" - qr) + (x' - q) 



+ ^ [Pr(a;" - qr) + P (g' - a;') + Sh] ■ 



(4.57) 



The derivatives are: 

_ 1 r/1-7 

(9g' 62 + ^ 



+ i 



<9pr/ „ X dqr 



, OSh , dSudq, 
+P + + 



(x - qr 



dq' dqr dq' 
1-7 



nigq + zrripg + (x -q 



1 + 7 

(x" - g,.) + (x' - g') 



1 
62 
1 

62 l\mqq + irriqp^ 

where we have used ^y- = — p', = Pr, and eqs. (|2.69| ). Similarly we find 



(4.58) 



dq' 

h 



1 



rrigg + im 



qp 



(x" - g^) - (x' - g') 



(4.59) 



The stationary conditions are satisfied if g' = x' and qr{q' ,p' , t) = x". This last equation 
defines p' implicitly, so that the contributing trajectory is the one that leaves x' at time 
zero and reaches x" at time t. The value of ^ computed at the stationary trajectory is 
simply iSn / fi. In order to perform the integrals we need the second order derivatives of 
^. The algebra is straightforward and the results are: 



1 2mqq + im 



qp 



62 TUqq + im. 



(4.60) 



qp 



iqp — 



m, 



qp 



pq 



h TTigq + im. 



(4.61) 



qp 
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= -4 • (4-62) 



Notice that all derivatives of ^ were calculated keeping the elements of the monodromy 
matrix fixed. Their variations would involve the computation of third or higher order 
derivatives of S. 

Inserting all these expressions into eq. ( |4.56| ), we find that the coordinate propagator 
becomes 



dgdP 



"IP 



^ 1 „iSH/fi 

^-^b~' ^= X _ (4.63) 

^jmqq + imqp I^^J ^qq^pp <,qp 



where Q and P are the variations of q' and p' from the stationary point. The quantity 
under the second square root in the last line is the determinant of the quadratic form 
in the second line. By eqs. (|4.60| ) to ( [4.62| ) it is equal to 

2 2 T"f\i QP (a \ 

KqqKpp ~ Knn ~ ToT I '■ N (^4.D4j 

''^ h\mqq + iniqp) 



and the final result is 



K^e,Ux\ t- x', 0) = ^ e^^«/^ (4.65) 



"qp 



which is Van Vleck's famous expression. 

Given the close similarity between Heller's wave-packet and ours, an identical 
calculation with ours will obviously give the result 

K{x", t; x', 0) = ^ ^ ^^{Sn+i)/n (4 56) 

h^j2Txi'mqp 

all classical quantities being calculated with the smoothed hamiltonian Ti instead of the 
Weyl H . This reduces to the Van Vleck formula for small and it might well be better 
for larger ^, but we have no evidence for this at the moment. 
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5. Comparison with the Herman-Kluk Propagator 



In addition to Heller's thawed gaussian approximation, there is in the literature a 
different initial value representation formula for the mixed propagator {x\K{t)\z) . This 
is the Herman-Kluk formula |Per84 |, derived by convoluting the Van Vleck propagator 
with coherent states and performing the resulting integrals by the method of stationary 
phase. This formula has been used many times in the last few years. 

We already pointed-out in section 4.1 the mistake of Grossmann and Xavier 
Gro98a|| in their tentative derivation of the HK formula from the semiclassical coherent- 
state propagator ( |2.59| ). In fact, in their original paper ||Her84|| , Herman and Kluk make 



a similar mistake in their evaluation of the stationary phase integrals. When performing 
these integrals they find (as we found in section 4.1) that the stationary trajectory is 
complex and given by mixed boundary conditions (equations (16) and (17) of ||Her84|| ) . 
However, instead of expanding the exponent of the integrand for this complex trajectory 
about a nearby real trajectory, they make a change of variables to initial position and 
momentum and assume these new variables to be real. Therefore, from the point of 
view of semiclassical analysis, the HK formula is incorrect. In subsection 5.2 we shall 
discuss why, in spite of that, it may still work. 

Here is the Herman-Kluk formula, for the mixed representation introduced in ([4.3| ), 
taken from ||Klu86|| with some adjustment of notations. 



{x\K{t)\z')HK = \/ ^ {nipp + nigq - inigp + tm. 



pq) 



exp 



(5.1) 



Once again, the term is absent and the classical hamiltonian used by HK is the Weyl 
hamiltonian. Besides these, there are two other differences between eq. ( [4.29|) and eq. 
( ^.1|) , both of them quite important at first sight. One is in the coefficient of the gaussian 
exponent, the other is in the prefactor. The coefficient (1 — 7)/(l + 7) of the gaussian 
in ( [4. 291) was essential for obtaining the right semiclassical widths ( [4.31| ) and ([4.32|) , but 
it is absent in HK! The difference in the pref actors is also astonishing. One formula 
has the square root occurring in the denominator and the other in the numerator. But 
actually the two prefactors are related by 



1 



[mpp + nigg - inigp + impq) 




(5.2) 



which shows that both differences between the two formulae imply that HK set 7 equal 
to 0. This is quite generally incorrect, except for the plain harmonic oscillator, for which 
both formulae give identical results. 

There is another major difference between formulae ( [4. 29] ) and ( p.l[ ), the 
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normalization, which is for (|5.1|) 

1 / f f (X - 

NuKit) = ^-^y^Y ("^gg + ^PpY + i^pq - ^qpY J dxexp I 

= , ^ • (5.3) 

This normahzation is not conserved since 7 is usually time-dependent. This is a grave 
flaw which shows up in almost every example. For the free particle, for instance, ( |4.29| ) 
is exact, while ( |5.1| ) fails completely in describing the spreading of the wave-packet. 

It may be argued with some validity that, since the |2;')'s form an over complete set, 
the fact that the propagation of each \z') according to HK is wrong does not necessarily 
mean that a wrong result will always occur when a |z')-integral is performed, as in eq. 
( [4.3|) . This is in fact a crucial point concerning the applicability of the HK formula, and 
we shall discuss it in detail in subsection 5.2. 

5.1. A numerical example 

We show here a detailed numerical illustration of the differences between HK and our 
semiclassical IVR formula. We consider the scattering of a particle with initial wave 
function {x\z') by a potential barrier. We choose the following Hamiltonian to test our 
results: 

H = ?^ + Vo [e"("-^) + e-"("+^)] 

2 (5.4) 

= |- + 2 K,e~"^coshaa; 

where Vq, a and A are parameters. The flrst term in the potential function represents 
an exponential wall located at x = +A. The second exponential closes the system at 
the left end with a second wall at x = —A. 

The wave-packet will be launched from a; = with positive momentum. If A is 
large compared to the coherent state width b, the packet will not see this left wall in its 
way towards the flrst collision. However, the fact that the motion is bound simplifles the 
quantum mechanical treatment, avoiding the complications of the continuous spectrum. 
We have set A = 5 and a = Vq = 1 for the potential and q' = 0, p' = 1 and b = 0.3 
for the initial coherent state \z'). Planck's constant was set to h = 0.05, which gives 
c = h/b^ 0.167. 

For short times the particle experiences almost no force, being well described by 
the free particle propagation. In this approximation the classical trajectory is just 

q = q' + p't 

(5.5) 

p = p 
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All the quantities entering the semiclassical formulae can be computed immediately and 
the result is 



'^qg '^pp 1 

Sh = p'^ t/2 - cH/A 
7 = itc/{itc + 2b) 

Substituting these expressions in [4.29| gives 



nipg = rrigp = ct/b 
Xr = cH/A 



(5.6) 



{x\K{t)\z') 



TT 40 2 

a/1 + itc/b 



exp 



1 1 / X — q' — p't^ 

2 1 + itc/b V b 



+ ^Ip' {x -q' - p't) + ]^p'q' + 



(5.7) 



which coincides with the exact quantum mechanical result. The HK formula, on the 
other hand gives 



{x\K{t)\z')HK = n'^b'^y^l -itc/2b exp 



1 X — q — pt 



+ ^lp'{x-q' - p't) + ^p'q' + p'h/2 



(5.8) 



This shows that HK not only describes poorly the width of the evolved packet but it 
also produces a non-physical increase in the total probability. 

For times of the order of A/p' the free particle approximation is no longer valid 
and we have to solve the problem numerically. The exact quantum mechanical solution 
of Hl'^n) = Enl'^n) for the Hamiltonian ( |5.4|) was performed by diagonalizing H using 
as basis states {x\(f)n) = sin {mix / 2L + n7i/2) where the width L is chosen such 
that V{L) = Emax- This is in fact the most accurate method (when you can use it). 
Emax is an upper limit to the eigen-energies to be calculated and it is related to the 
number of basis states used in the diagonalization by E^ax = A^^vr^/i^/(8L^). This 
guarantees that the basis states span the energy interval from to Emax and that the 
part of the potential with —L < x < L has V{x) < E^ax- We have used = 400 in 
our computations, which gives Emax ~ 9.4 and L 7.2. With this choice the first 260 
energy levels converge with at least 5 digits (as compared to a larger diagonalization), 
spanning an energy interval from to 5.5. The initial wave packet {x\z') can be easily 
expressed in terms of the basis and, therefore, in the basis |\E'„) of eigenstates of H. 

The semiclassical calculation of both HK and our semiclassical formula needs only 
one single trajectory starting from {q',p') and its tangent matrix to compute the whole 
function {x\K{t)\z'). In our formula this trajectory is a solution of Hamilton's equations 
for the smoothed Hamiltonian {v\H\u), which is given by 



2 Vie""^ coshax + h'b'/A 



2 1,2 



(5.9) 
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where Vi = Vq exp (6^a^/4). The resuhs are shown in Figures la to Id. Figure la shows 
the square modulus of the wave function for the times t = and t = A. The gaussian at 
t=0 represents the initial coherent state. The dashed line shows the (properly scaled) 
potential function and the dashed vertical bar shows the location of the classical turning 
point. The lines at t = 4 show the time evolved wave-function according to the exact 
calculation (solid), our formula (dashed) and the HK approximation (dotted). Notice 
the increase in the height of the HK gaussian, in opposition to the exact evolution, 
which spreads and decreases. Our semiclassical formula is in very good agreement with 
the exact calculation. For t = 6, Figure lb, one has the impression that HK is not so 
bad, since the wave packet height increases again. That is, however, only a fortuitous 
occurrence, as demonstrated in Figure Ic, for t = 8 and Figure Id for t = 10. For longer 
times the peak of the HK gaussian increases more and more, whereas the exact wave 
packet height decreases to compensate for the spreading. At times of the order of 25 
(not shown) the exact propagator (and our semiclassical formula) has a peak of height 
around 0.2 while HK's peak is around 8. At this time the width in HK is also completely 
wrong. 



5.2. Van Vleck's formula and the HK propagator 



We saw in the previous subsection that the propagation of gaussian states via Herman- 
Kluk is very imperfect. It works only for harmonic potentials or for very short 
propagation times. In spite of this, the HK formula became rather popular in the last 
five years, and was applied to several problems in Chemical Physics. Among those we 
cite the photodissociation of CO2 [ Wal95 |, the colinear scattering of H2 by H [ Uar9(: 1, 
and the non-adiabatic dynamics in pyrazine molecules ||ThoOO|| . Several theoretical 



articles were also published, testing the HK formula in model systems |[Kay94b| , |MaiO0| 



or proposing modifications in the formula for increasing its accuracy ||Her97| , p?hoOO|] . 
All these applications deal, not with the plain mixed propagator {x\K(t)\z') , but with 
integrals over the phase space variables q' and p', such as in eq. (^l3| ). 

The question is then how a poor time evolution of gaussian wave-packets, as in 
HK, may produce acceptable results after integration over q' and p'. In this subsection 
we shall clarify this point and, at the same time, show why HK is so expensive 
computationally ||Wal95|| and sometimes just does not work [pVIccOOj] . First of all we shall 



calculate the coordinate propagator ([4.55|) using the semiclassical formula of Herman- 
Kluk for {x"\K{t)\z'). Once again the integration over q' and p' will be performed by 
the stationary exponent approximation, just as we did in subsection 4.4. We shall see 
that HK does recover Van Vleck's formula ||Van28|| exactly! And that is the reason why 
sometimes it works. At the end of this subsection we shall comment on this result and 
discuss why HK does not lead to long-time convergence when the integrals over q' and 
p' are performed numerically. We shall also point out why the HK prefactor diverges 
for chaotic trajectories when the correct prefactor should, instead, go to zero. 

We use the same notation as in subsection (4.4) and follow the same steps. The 
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coordinate propagator for HK can be written 



KHK{x",t;x',0) = J -|^^ ^ y ^ {rupp + rriqq - irugp + irripq) S'^'^^"^ (5.10) 
where the exponent ^{q',p') and its derivatives are: 

^ = - ^'■)' + - ^')'] + - ^'-) - + (5-11) 

= ^[(^" - ^-)K'? + ^"^P^) + - ^')] (5.12) 
= -^[{x" - qr){mpp - irrigp) - {x - q')] . (5.13) 

The stationary conditions are satisfied if q' = x' and qr{q',p', t) = x" and the value of ^ 
computed at this trajectory is iSu/h. The second order derivatives are 

i 1 
iqp = ^pq = --^[1 - rnqqijUpp - irriqp)] = --mqpirriqq + impq) (5.15) 

ipp = -^rnqpirupp - imqp) . (5.16) 

Inserting all these expressions into eq. (|5.10|) , we find that the coordinate propagator 
becomes 



KHK{x",t;x',0) = 71 h ^\h:{mpp + 



2 

dQdP 
— exp 



li^qqQ'+2^qpQP + ^ppP') 



iSu/h 



-1,-1 ^ ^'^"'^ 

TT 26 . {^rripp + mqq - luigp + irripq) x == (5.17) 

* l^ykqq^PP ^qp 

Once again the quantity under the square root in the denominator is the determinant 
of the quadratic form in the second line, which is 

^qq^pp ~ ^qp — '^f^qpij^pp + f^qq ~ ilTT'qp + ilTT'pq) ■ (5.18) 

The final result is 

Khk{x", t- x\ 0) = ^ e*^«/^ (5.19) 

liiimqp 

and it coincides with the Van Vleck formula ( |4.65| ). 

We can now understand why HK may sometimes work, even if it is not a correct 
semiclassical formula, and we can also understand the origin of its main drawbacks. 
The latter include the very slow convergence of the integrals over g' and when 
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done numerically |[WanOC , MccOO |, the lack of normalization | Her84 1 , and the blow- 
up of the prefactor for chaotic trajectories ||MaiOO|| . In fact, numerical convergence 
of the integrals over q' and p' is often achieved only after resorting to smoothings 
Her 9 7] , [MccO0| , fThoOOl , |SunOO|| , and the results have not always been satisfactory. 



The reason why HK does sometimes give good results is precisely that it is able 
to recover Van Vleck. In fact, Kay ||Kay94a|| derived the HK formula by making an 
ansatz for the propagator and demanding that this ansatz satisfy the basic condition of 
agreeing with Van Vleck. In words, if one insists in keeping the width of the propagated 
wave packet constant, then one must arrive at the HK pre-factor if one wants to get Van 
Vleck when integrating over q' and p'. The price paid for doing so is non-conservation of 
the norm, a high price to say the least. Therefore, if HK is to be used at all, it has to be 
under an integral. However, even when integrated, HK leads to very slow convergence 
and oscillatory behavior, especially at long times. To see why this happens, we note 
that the square root of the determinant of the quadratic form in the stationary integral 
( p. 171 ) is a measure of the phase-space area Aq'Ap' that matters in the integration. We 
have 



{Aq'Ap')HK 



1 



n 



1 



m, 



qp\ 



(5.20) 



As we saw in subsec. 5.1, \J\mpp + m^^ — iniqp + impq\ increases with time, for our 
numerical example as well as for the free particle. This means that the relevant phase- 
space zone of initial trajectories that one needs to sample, gets smaller with time. The 
physical interpretation of this is very simple: any initial swarm of trajectories spreads as 
time passes. However, only those trajectories in the neighborhood of the one connecting 
x' to x" contribute significantly. The size of this neighborhood is determined by the size 
of the propagated wavepacket. In the case of HK the propagated packet keeps its width 
fixed at AqAp = be = h. For long times the initial spread of trajectories that end inside 
this small region shrinks very fast. This is why the numerical integration of eq. ( |5.1CI| ) 
by sampling initial trajectories is bound not to converge for long times, since very few 
trajectories are going to be picked up in the relevant region. 

Compare this with the relevant phase-space spread for the integration with our 
formula 



{Aq'Ap' 



1 



h 



iqqipp if 



qpi 



2m 



qp 



miqq -\- tmiqp | 



(5.21) 



We know that 



m. 



qq 



im 



qp\ 



increases with time, since the wave packet spreads, and so 



does the phase space region of contributing trajectories, ensuring an integration which 
must be more efficient. 

We also mention that the application of semiclassical formulae to chaotic systems 
is known to be difficult due to the exponential proliferation of contributing trajectories 
from x' to x" for long times. This is compensated in part by the exponentially small 
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contribution of each individual orbit to the propagator. The Herman-Kluk prefactor, 
however, assigns an already divergent contribution to each of these trajectories, leaving 
no hope for accurate results [[MaiOO| . 

As a final comment we note that the smoothing technique introduced by Herman 
Her97|| has the role of cutting off the contributions from trajectories whose action 



have a large first order variation. In our semiclassical formula this cutoff is performed 
automatically by the prefactor. This is in fact exactly what one expects from a stationary 
phase integral. 
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Figure 1. Square modulus of the time evolved coherent state launched from q'=0 
and p'=l. The dashed line shows the potential function and the dashed bar shows 
the location of the classical turning point. The solid, line shows the exact propagation, 
the dashed line shows our semiclassical formula and the dotted line shows the HK 
approximation, (a) initial wave packet att = and exact and semiclassical calculations 
att = A. (b) exact and semiclassical calculaiions att = 6; (c) exact and semiclassical 
calculations att = 8; (d) exact and semiclassical calculations at t = 10. 

1 T7^Kv.n,„v.Tr onnc 



The energy representation 



52 



6. The energy representation 

In the special case of time-independent systems, the Fourier transform in time of 
the coherent state propagator, which we shall refer to as the (coherent state) Green's 
function, has poles at the quantized energy levels. For its diagonal elements {z' = z") 
the residues are the so-called Husimi distributions, i.e. the absolute squares of the 
eigenfunctions in the coherent state representation, which are also called the Bargmann 
wave functions. In this section we shall derive the semiclassical expressions for the 
poles and residues of the Green's function. We shall obtain from them the semiclassical 
quantization rule for the energy levels and the semiclassical Husimi distributions. 



6.1. The monodromy matrix in one dimension 

When the trajectory is a periodic orbit, with the initial point returning to itself after 
the period T, the tangent matrix M of section 2.6 is called the monodromy matrix. In 
the next subsections it will be necessary to know the element for a real periodic 
orbit traversed n times, i.e. for t = uT . We shall calculate it here. First we shall 
review a general formula for arbitrary canonical variables, then we shall switch back to 
the (m, v) variables. 

Let Q and P be canonical variables and H{Q, P) the Hamiltonian. We want to 
find the monodromy matrix M for periodic orbits in the {Q, P) representation. M has 
4 elements, which we shall determine from 4 linear equations. For the first 2 equations 
we apply M to {Q, P), the velocity vector in phase space. This should yield this vector 
again, because a small displacement along the trajectory maps onto itself. Thus 

For the other 2 equations we consider a small displacement {6Q, 6P) non-colinear with 
[Q, P), and we assume for now that the Hamiltonian is not harmonic, so that the period 
T of a periodic orbit depends on its energy. Then {6Q, 6P) points to another periodic 
orbit, with a different period T + 6T and a different energy E + 6E. When we apply 
the M matrix to {SQ, SP), since the propagation time is only T, the result falls short of 
the original displacement by 5Tx (velocity), hence 

-5T{^] . (6.2) 





Obviously, given the 4 equations ( |6.1|) and (6^), M must have the form 



M = l + N (6.3) 
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where N must obey the 4 equations 

^(11=0 (6.4) 




N Tr. =-ST I . (6.5) 



From eqs. (p.4|) one sees that must have the form 



Pb -Qb 

with a and b still unknown. From eqs. (|6.5| ) one calculates a and 6 as 



where Ji?, the energy difference, actually appears as Q6P — PSQ, which is 6E 
by Hamilton's equations. The final form of the monodromy matrix in the {Q, P) 
representation is therefore 



(6.8) 

\P' -QP) 

where dT/di? is the derivative of T{E), the period of the orbit as a function of its 
energy. Although we restricted ourselves to the anharmonic case, it is obvious that this 
formula holds for the harmonic case also, for which dT/dE vanishes. 

We shall now apply this to the variables u and v. But the latter are not quite 
true canonical variables because of the factor ih in the equations of motion ( |2.35D . To 
avoid this source of confusion we define temporary canonical variables U = y/ihu and 
V = y/ihv. Then the general result ( |6.8|) certainly applies to U and V. When we 
reintroduce u and v, we find for the monodromy matrix in the {u,v) representation 

M{u, v,T) = l+ ih—r (6.9) 
dE 

where F is given by 

r=h . (6.10) 

It is easy to check that, because F^ = 0, the monodromy matrix for n traversals of the 
orbit is just 

irri 

M{u, V, riT) = M{u, v, T)" = 1 + nih—T. (6.11) 

dE 

In the next subsections we need only one matrix element, which we denote by 

(M")^^ = 1 - nih—uv. (6.12) 
dE 
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6.2. Prologue: The Green's Function for the Harmonic Oscillator 

The Feynman-type propagator Kf{z" ,t; z' ,0) in the coherent state representation, 
hke other such propagators in quantum mechanics, is discontinuous at t = 0. One 
distinguishes between the "forward propagator" , which is defined by eq. ( p.5| ) for t > 
but vanishes for t < 0, and the "backward propagator" which vanishes for t > 0. The 
Green's function is the Fourier transform of the forward propagator, hence it is obtained 
by integrating the time only in the interval t G [0, oo]. We are interested in its diagonal 
elements, which we denote by 

G(z, E + i-f):=— / AtKiz, t- z, o)e^(^+*^)*/^ . (6.13) 
in Jo 

The small positive quantity 7 is introduced for convergence; at the end it will be set 
equal to 0. When we insert the semiclassical propagator ( p.68|) , G becomes 



1 1 
G{z, E + t^) = -J2 dt^=_ exp IMt)] (6.14) 
V Jo V(M,)„„ 



with 

Mt) = I 



S^iv", u\ t) + X,{v\ u\ t) + j (l^'f + \z'\^) + {E + i^)t 



(6.15) 



Instead of the variables v" and m', we could also have written eq. ( |6.15D in terms 



of the single variable z, since u' = z and v" = z*. But this does not mean that the 
trajectory going from z to z in some arbitrary time t is periodic. As we have already 
discussed (see before eq. ( p.37|) ), the trajectory is usually complex and neither u nor v 



match at the end points. However, for the case of interest here, namely one-dimensional 
bound systems, there exists for each phase-space point z a (minimal) time T = T{z) for 
which the orbit through z is periodic and real. Then the end points do match and T is 
the period. Obviously this particular time t = T will play a special role in the integration 
of eq. ( |6.13|) , and so will its multiples t = nT, which correspond to repeated traversals 
of the periodic orbit. Therefore, after evaluating the integral ( |6.13| ) by the stationary 
exponent method, we shall expand the stationary time to = to{z, E) about the classical 
period nT{z). The resulting Green's function will be a good approximation to the actual 
Green's function G{z, E) only for arguments z and E satisfying to{z, E) ^ nT{z). 

Before we carry out the calculation of G{z, E) for a general Hamiltonian, it is 
instructive to find out the result of these approximations, stationary exponent plus 
expansion about to = nT, for the harmonic oscillator. This will help us understand the 
nature of the approximations. 

We consider an oscillator of unit mass and frequency uj, and we take the parameter 
b of eq. ( |2.3| ) to be ^Jthjuj. Then the semiclassical propagator is easily shown to be 

K{z, t) = e-'"*/2 exp \\zWe-''^' - 1)1 , (6.16) 
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which is also the exact result. The first exponential is the prefactor M^^^"^ . This K is 
now carried into eq. (|6.13|) , which becomes 

G(z, E + ij) = - / dte-*-*/2g.(£;+n)V^H-kP(e~'--i) . (6.17) 
^ri Jo 

The whole of the second exponent is proportional to 1/h, when is written in terms 
of q and p. But there is no h in the prefactor. Hence we are exactly in the situation 
mentioned at the end of section 3.3. We look for the time to at which the second 
exponent is stationary. We find that to is given by the condition 

E + i7 = ^e-^'^*«|z|2 . (6.18) 

The solutions of this equation can be written to = + "^T, where n is an integer, 
T = 2tt/uj, and Tq = To{z,E) is pure imaginary (in the limit 7 0). Thus there is 
an infinite number of times at which the second exponent is stationary. For each to 
we expand this exponent to second order in the vicinity and we perform the gaussian 
integral. Then we add all these results together to obtain the following G 

G(Z E + i-f) = — Y^ g-f (nT+To)gf(i?+i7)(nT+To)+|^|2(e--^0-l) / ^_ ^ (t_to)2e--To 

'^tl Jo 

= ^ / ^^(E+i^-!uo/2) g|^|2(e— l'o-l)+i(E+»7-;K^/2)To/ft 

ih\z\\J u;2g-ia;To 



n=l 



1 I 27T etf(^+^7-?i-/2) 

e 



\z\^{e-^^'^0 -l)+i{E+i'y-hw/2)To /H 



(6.19) 

In the above we started the sum over n at n = 1. Why not n = 07 The answer is 
that it does not make any difference, as long as we use this theory only to calculate the 
energies of the stationary states and their Husimi distributions. This is discussed in the 
paragraph following eq. (|6.34| ) in the next subsection. 

The poles of ( |6.19| ) can now be found in the limit 7 — 0. They are given by the 
condition E =: E^ = hjjj{m + 1/2) with integer m. The residue for E = E^, which is 
the Husimi distribution, is 

PnusU^) = , I el^l^(-'^^-^)« . (6.20) 

Using (|6A8| ) for To with E = E^ and 7 = 0, we obtain e-'^'^^\z\'^ = m + 1/2 and (lOOp 
becomes 

1 / Ul2 

^m+l/2-|2|2 



v/27r(m + 1/2) \m + l/2j 

This result should be compared with the exact Husimi distribution 

_ 1 -|.|2|^|2m ^ (6.22) 
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which is an annulus with its maximum at = m. The mean value of is m + 1 
and the logarithmic second derivative at the maximum is — 1/m, giving a width of order 
^/Tn. Both formulas, exact and semiclassical, have the same 2;-dependence, namely, 
|2;|2"ig-|2p_ They differ only in the coefficient. Obviously ( |6.22 ) is correctly normalized 



and ( |6.21| ) is not. The difference is small, however. Using Stirling's approximation for 
ml, we can write ( |6.21|) as 



-1 / 1 \ '"+1/2 

PHusimi(^) = ^e-l^lVP'" 1 + T^—r . (6.23) 

ml \ 2m + 1 / 

For large m the quantity to the left of e~^/^ becomes e^^/^. Thus the semiclassical Husimi 
becomes exact in the limit of large quantum numbers, and it is a valid approximation 
for all phase space points z. 

For general Hamiltonians we shall not be able to do such a complete calculation. 
As we shall see, the stationary time to is given by an implicit equation that cannot 
be generally solved. The best thing to do to get an explicit result will be to expand 
to about the classical period nT. Now we can check this additional approximation 
explicitly for the harmonic oscillator. Let us pretend that eq. ( |6.18|) cannot be solved 



exactly for tQ{E, z) and, instead, let us solve it by expanding to about nT. So we write 
to = nT + To and consider only small values of To. Taking 7 = in eq. (|6.18| ) and 
expanding for small Tq, we find 

(6.24) 



hjjj 


z 


2 




E 




z 


2 



Since To is independent of n we may go through the same steps as in eq. ( |6.19| ). The 



position of the poles does not depend on Tq and it is not affected by this approximation. 
The Husimi distribution, however, does depend on To explicitly. For E = E^ we get 

(6.25) 



\z 


2 

— m 


-1/2 




z 


2 



To is small whenever |zp is close to m + 1/2, which is the classical orbit with quantized 
energy Em- Therefore, for each eigenstate, there is a phase space region centered on this 
classical orbit where the approximation is justified. This region should encompass most 
of the distribution when m is sufficiently large. Expanding eq. ( |6.20| ) to second order in 
To and using ( |6.25| ) we obtain 



1 r bp -m- 1/2 2] 
PHusimi(2;) ~ , — , , exp — — . (6.26) 

The maximum of this annulus comes aX \z\^ = m + 1/ 8m + ■ ■ ■ , which agrees with the 
exact maximum for the first two orders of 1/m. The logarithmic second derivative also 
agrees for the leading order. Therefore, the result of expanding the stationary time 
about the classical period is to restrict the region of validity of the Husimi distribution 
to not too small quantum numbers. 
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6. 3. The Green's Function for General Hamiltonians 



We now return to our general calculation, eq. (|6.13| ). In what follows we shall omit 
the subscript z/. It will be replaced shortly by the multiple traversals around periodic 
trajectories. The stationary exponent condition is given by 0'(to) = 0, with given 
by ( I05D . It is 



0'(to) = 

The third equation ( p.56|) gives 




to 



dl 



+ E + ij 



0. 



(6.27) 



to 



dS/dt\ 



t=to 



-£{z,to) ■ 



(6.28) 



This script £ is the energy of the classical trajectory, not to be confused with E, which 
is the energy variable in the Green's function. For the first time derivative of X, defined 
in eq. ( p.61|) , we have 



dl/dt\ 



1 d^n 



t=to 



2 dudv 



t=to 



1 

+ 2 



to 



d^n dz 



d'^zdz* dt 



+ 



d^n dz^ 



t=to 



dzd'^z* dt 



t=to 



dt' . (6.29) 



The terms under the integral involve third derivatives of Ti. with respect to the phase 
space variables. We have consistently neglected such terms up to now and we must do 
so here again. Hence we use only the first term, to which we give a simpler name 



dl/dt\,__ 



1 d^n 



=*« 2 dzdz* 
The stationary condition is then written 

E + - £{z,to) + e{z,to 



e{z,to) 



. 



(6.30) 



(6.31) 



Besides the contributions from the stationary points t = to, there is also a contribution 
to G coming from the vicinity of t = 0. As we shall see later, this contribution is not 
needed in the calculation of the energy levels and the Husimi distributions, hence we 
shall ignore it. We have also discarded it in the previous calculation for the harmonic 
oscillator. 

In order to perform the time integral, we need to expand 0(t) around to to second 
order. The first derivative at to is zero by definition. The second derivative is 



0"(to) 




+ 



to 



dt^ 



(6.32) 



to. 



For the second derivative of S we introduce the notation 

a{z,to) := d^S/dt"^ 



\t=to 



(6.33) 
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From ( 6.29|) we see that d'^l/dt'^ involves only third or higher order derivatives of Ti 



and therefore we can discard it completely. 

We can now calculate G{z,E + ^7). The effect of our expansion is that we have 
applied once again the gaussian approximation to the integral in eq. (|6.13|) , which is 
then straightforward: 



G(Z F,^i^) = ± ^ ^i[S{z,to)+Iiz,to)+m^\' + iE+i^)to] ^f%M(t_t„)2 



27tih ( i 



ihy M^^{to)a{z,tQ) 



exp I ^ [S{z, to) + I{z, to) + z^lzp + {E + i-f)to] | 

(6.34) 



Actually this integral is a sum of gaussians, because there are many solutions to eq. 
( |b.31| ), as we have seen with the harmonic oscillator. As we mention later, a has a 
positive imaginary part for to ~ nT; this follows from eq. (|6.53|) , plus the fact that S 
is real. It means that the integral is strongly convergent for very large times, both 
positive and negative. However, the integral in eq. (|6.13| ) does not extend over the 



interval —00 < t < +00, but over < t < +00: is there a problem at t = 0? Yes, 
there might very well be a problem, both for t = and for to ~ nT and small n. 
There is no problem for large n, since nT becomes arbitrarily large. We intend to use 
this Green's function only for the semiclassical calculation of the energy levels and the 
Husimi distributions. These are determined by the poles of the Green's function and 
by their residues, respectively. And the poles and residues are determined solely by the 
behavior of the integrand of ( |6.13| ) at very large times. Hence we are safe in ignoring 
possible mistakes at small t. In particular, the contributions from trajectories with 
t ^ give rise to the Thomas-Fermi approximation, which we shall not consider here. 
Therefore it is all right for our purpose to calculate the integral as if it extended from 
—00 to +00. 

We shall now transform the prefactor by finding a convenient way of expressing 
a{z,to)- In the many differentiations which follow, it is important to remember that 
the independent variables in 5* are u',v", and t . One should do the differentiations 
first, and only afterwards may one compute the functions for the stationary orbit by 
replacing t hj to, u' by z and v" by z*. First we write two different ways of expressing 
dS/dt 

dS/dt = -n [u', v'{u', v", t)] dS/dt = -n [u"iu', v", t),v")] . (6.35) 
This leads to two different ways of writing a 

a{z,to) 



a{z,to) 







dndv' 


dt^ 


t=to 


dv' dt 






on du" 


dt^ 


t=to 


du" dt 



t=to 



(6.36) 
(6.37) 



t=to 
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In eq. ( |6.36|) we use Hamilton's eqs. (|2.35|) to express dH/dv' and the second eq. (|2.56| ) 
to express dv' /dt . In eq. (|6.37|) we do similar transformations with the other variables. 
This gives the two forms 



ihu' 



h du'dt 



ihv" 



t=to 



h dv"dt 



(6.38) 



t=to 



Referring now to eq. ( |6.35| ), using Hamilton's eqs. once again, and using the first and 
second eqs.( |2.56| ) once again, we transform the second derivatives of S as follows 

d^S dniu",v")du" 



du'dt 



du" du' 
dn{u',v') dv' 



— %fiD^^ 

h du'dv" 

-ihu'—- 



(6.39) 
(6.40) 

dv"dt dv' dv" h du'dv" 

After cancellation of the Ws, both forms of the equations give the same result for a, 
namely 

d^S 



a{z,to) 



du'dv" 



u'v" A, 



t=to 



\t=to 



ihu'v"/My^{z,to) 



(6.41) 



Amazingly, the denominator under the square root in eq. ( |6.34| ) has become simply 
ihu'v", where the velocities are computed at the complex trajectory with t = to- One 
should, however, be careful about the phase of the M^y under the square- root, as already 
mentioned in sections 2 and 4. We shall discuss it in a moment. 

In spite of this great simplification it is still hard to find the poles of G{z, E) by 
looking at eq. ( |6.34D . This is because the stationary time to is given implicitly by eq. 
( |6.31| ) and it refers to complex trajectories. According to our discussion following eq. 
( |6.15| ), large contributions to G{z, E + i'j) are expected for to close to nT, where T is the 
period of the real orbit through z. Therefore, as we did for the harmonic oscillator, we 
proceed to expand eq. (|6.34|) about to = nT, summing over n. We write to = nT + To, 
substitute it into eq. ( |6.31| ), and expand the terms to find To: 



E + i-f- £{z, nT) + e{z, nT) - d£{z, t)/dt\^^^^ To + dt{z, t)/dt\ 



t=nT ^0 



. (6.42) 



Equation ( |6.30| ) shows again that de{z,t)/dt involves only third or higher derivatives of 
Ti. and therefore we discard it. For real periodic orbits, the functions S and e depend 
only on z and we have 



d£iz,t)/dt\^^^T = - d^S/dt' 
Equation (|6.42|) can then be written 



\t=nT 



-a{z, nT) 



-a'-'^Hz) 



{E + i-f - S{z,nT) + e{z,nT)) +a^'''>{z)To ^ . 



(6.43) 



(6.44) 



Our basic assumption is that Tq ^ nT for each solution to labeled by n. Hence the sum 
of the terms inside the parenthesis above should be small. Solving eq. ( |6.44 ) for Td gives 

E + i-f-£{z)+eiz) 



To 



a 



(6.45) 
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In the limit 7—^0, the points z satisfying £{z) — e{z) = E have to{z,E) = nT{z). 
When we expand the Green's function about to = ^T, we are restricting the vahdity 
to the neighborhood of these points. Notice that, contrary to the case of the harmonic 
oscillator, Tq does depend on n and the sum over multiple traversals has to be performed 
after the expansion about to = nT. 

We are now in a position to discuss the pre-factor which, according to eq. ( p.41| ), is 
given (up to its phase) by the square root of u'(to) ''^"(^o)- Its main contribution comes 
from periodic orbits and, once again, we have to expand each of these velocities for to 
close to riT. Using Hamilton's equation for 'u'(to) we get 



M'(tc 



h dv' 



t=to 



To 



To 

t=nT 



t=nT 



h dv'"^ \ h dtdu' 
z + 'jH,, z* To/ Ml 

z + ^n,. To (6.46) 



where we have used eqs. ( ^.56| , |6.39| , ^.67] ) and ( |6.41[ ). Doing a similar calculation for 
v"{to) we find 

u'ito) v"{to) ^ (1 + 0(a(")To)) . (6.47) 

Now we want to argue that the terms of order a^^^^To are small and that, in fact, the 
prefactor can be simply calculated at the periodic orbit itself, so that ( |6.47| ) is just 
lip. We gave a similar argument in section 3.3, where the neglected terms were of 
order h compared to the terms kept. Here, however, their order is a/^, not as small, 
but still small! We shall verify this a posteriori once we have calculated the Husimi 
distribution. We shall find later that the Husimi attains significant magnitudes only 
when E — £{z) + e(z), which is the same as a'^'^^To by ( |6.45| ), is of order \/h. 



To calculate the phase of the prefactor, or the tangent matrix, for a real periodic 
orbit, one must follow it when moving around the orbit. The harmonic oscillator provides 
again a simple illustration of what happens: in this case the solution of Hamilton's 
equations ( p. 351) with initial conditions m(0) = u' , f (0) = v' is simply u{t') = u'e^^^^' 
and f (f) = f'e*'^*'. The tangent matrix is diagonal and My^ = e*^*. The pre-factor of the 
time-dependent propagator is, therefore, e"*"^*/^ (see eq. (|6.16|) ). When this is calculated 
at a periodic orbit, we see that after one period, for t = T = 27r/ci;, Myy has rotated by 
27T and the phase of the pre-factor is e~*'^ = —1. The phase for t = 2mx /uo is just e~*"'^. 

This phase of — nvr is not particular to the harmonic potential; it is a consequence 
of the fact that the motion is periodic in two-dimensional phase space. To see this, 
consider a periodic orbit of energy £ and period T in a generic system. The tangent 
matrix propagates small displacements about this orbit. Any such small displacement is 
going to point to a nearby periodic orbit with energy £ + 5£ and period T + 5T. As we 
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move once around the orbit the displacement vector stays hooked to these two reference 
orbits and, therefore, has to rotate once as well. The total rotation is not exactly 27i, 
since the period of the nearby orbit is slightly different from T. The angle of rotation is 
actually 271 + 6, where the 6 piece follows from eq. ( |6.12|) for M^y. When M^^ is raised to 
the —1/2 power in the prefactor, this 2tt contributes a phase of — vr to the propagator, 
just as in the harmonic oscillator. Then there is the angle 6, which is a function of n 
and which remains in the prefactor. Equation ( p.l2| ) shows that 6 never gets very large; 
for all n's it is bounded either by and 7r/2, or by —7c/2 and 0, depending on the signs. 
However, by virtue of eq. ( |6.41| ) the phase 6 disappears completely from the calculation 
and the only relevant phase from the square root is —nir. From now on we shall take 
this phase out of the prefactor and include it explicitly in the exponential. Replacing 
u'v" in eq. ( |6.41| ) by |ip we get 



G{z, E + ij) = exp I - [S{z, to) + I{z, to) + ih\z\'^ + {E + ij)to - nnh] | . 

(6.48) 

Next we expand 5" and I around to = riT. Since the zero-order orbit is the real 
periodic one, we may write, using eq. (|2.52|) 

S{z\ z, nT) = nS{8{z)) - nS{z) T - ih\z\'^ . (6.49) 

Here S{z) is the energy of the real periodic orbit going through z. It should not be 
confused with E, the quantal energy variable. And S{S{z)) is the "reduced action" 

ih f'^ 

S{£{z)) = — {uv- vu)dt' (6.50) 



2 Jo 

which, for the real periodic orbit, is a function only of S, or of the period T, themselves 
functions of z. The third eq. ( [^.56D gives 

dS/dt\,^^^ = dS/dt\,^^ = -S{z) . (6.51) 

Either £ oi T can be used as the independent variable for the periodic trajectory. The 
second derivative of S becomes (see eqs. ( |6.33| ) and ( |6.41| )) 

d''S/dt\^^^ = a^^\z) = 2^|i|V^i"^ • (6-52) 
Using eq. (|6.12| ) and dS /d£ = T, we find 



1 Mi"^ 1 n dT . . 1 d^^ 



a 



W ih\z\'^ ih\z\'^ \z\^d£ ih\z\^ d£ 



n—- . (6.53) 



From the definition ( |2.61| ), we have 2{z* , z,nT) = nl{z* , z,T). Once again, for the 
periodic orbit, X depends only on the classical energy or on the period and we can call 
it I{£{z)). For the first time derivative we have by ( |6.30D 



dl/dt\,^^^ ^ 1^ = e{z) . (6.54) 
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Calling ijj^E^z) the exponent in eq. (|6.48|) and using eqs. (|6.49|) -( |63^ we get, up to 
second order in Tq, 



/ 1 



S{8{z)) + X{8{z)) -nh+^{E- S{z) + ^7) 



{E-S{z) + e{z)+z^)n + ^T, 



(n) 



(6.55) 



Next we add and subtract e d5/d£^ inside the brackets on the first line and use eq. ( |6.45| ) 
for Tq in the second line. The term i 
When we add these two terms we get 



for Tq in the second line. The term in Tq becomes minus one half of the term in Tq. 



i^iE^z) =-n 



2;ia(") 

Finally we use eq. ( p.53| ) for a^"-' 



S{£{z)) + X{£{z)) -^^+%{E- £{z) + e{z) + .7) - 
" (E-8{z)+e{z)+^^)\ 



(6.56) 



'^{E,z) =-n 



S{£{z)) + l{£{z)) -^h+^{E- £{z) + e{z) + ^7) - e(^)§ 



I (PS 



2h^z\ 



■{E-£{z)+e{z)+i^f . (6.57) 



The first, fourth and sixth terms in ( 6.57 ) are the Taylor expansion of the function 
S{E + + e{z)) to second order around £. This is an acceptable approximation as long 
as E — £{z) + e{z) + 27, i.e. Tga'-^^ is small. Since e is of order h we can further write 
S{E + 27 + e{z)) ~ S{E + 27) + T{E + i'y)e{z), neglecting terms of order h"^. Then eq. 
( |6.57| ) can be written 

ij{E, z) ^ ^jn{S[E + «7) + Z[£) + e(z) \T{E + 27) - T{£)\ - ixh] 

Since X and e are themselves of order we may also replace £^ by ii^ + Z7 in the 
argument of both X{£^ and T(£\ the error in doing so being of order T^a^'^^fi. This 
simplifies ip{E, z) even more and we get just 

1 



^{E,z) 



h 



n 



{S{E + ^7) + X{E + ^7) - Tih] 



2h^ 



— {E-£iz) + e{z)+z^y. 

(6.59) 



Inserting this as the exponent in eq. 



and summing over n we obtain 



G{z,E + i-f) = -T^X^expi^n [5(^ + ^7) +X{E + t-f) - nh] 



n=l 



{E-£{z) + eiz)+t^y 



2n2|iP 



(6.60) 
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We may now perform the sum over all multiple traversals. The result is 



+ = i_e-(^+x-..)M ^^P| ^2 > (6-61) 

where S and X are taken at E + i'j. 



6.4- Semiclassical energy levels and Husimi functions 

Since the poles of the exact G are the stationary state energies, and their residues are the 
corresponding Husimi distributions, we are now in a position to calculate semiclassical 
approximations to these quantities. We let 7 — 0. If 7 7^ the poles are displaced by 
27. Poles occur whenever {S + X — irfij/h = 2m7r. This is the quantization rule, which 
can be rewritten 



{S + X){EJ = {m + l/2)h . (6.62) 

This is what replaces the usual WKB formula, which can be obtained from the 
coordinate Green's function. One should not forget that S here is not the same as 
the S in the WKB formula, and that the latter does not contain any X, of course. One 
should also realize that the +1/2 is obtained in very different ways in the two cases. In 
the usual WKB, the semiclassical approximation diverges at the turning points, and one 
must use some delicate arguments around this problem to derive a connection formula. 
Here, on the other hand, there is no divergence and we derived the +1/2 with a simple 
continuity argument. 



The residue of 



at = is 



i gi(5+X-7rft)/ft 
^ I _ Qi{S+I-nh)/h 



1 1 



diS + X)/d£ T(E„)+ (dJ/d^)|^,^ • 
Hence the residue of G, which is the Husimi distribution for the level Em , is 



(6.63) 



(6.64) 



j -{Em-S{z) + e{z)y ^ 

^^--^^^ = \z\ [nEm)+{6X/d£)U 1 ^ • ^'-''^ 



It is centered close to the classical trajectory with the quantized energy. Both its 
amplitude and its width are modulated by the phase space velocity \z\. The width is 
given by the gaussian factor. In terms of either variable £ or z, it is of order ^/h. This 
follows from the fact that is of order unity (see eqs. ( ^.2[ ) and ( ^.3] )). Therefore, 
our decision not to expand the pre-factor (see eq. ( |6.47D ), and the approximation after 
eq. (|6.58|) , amount to discarding corrections proportional to ^/h in the Green's function. 
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These are not as small as the other corrections to the semiclassical approximation, which 
were of order h (see section 3.3), but they are still small and vanish in the semiclassical 
limit. 

For the second type of path integral discussed in section 3, the quantization rule 
can be obtained by changing I to —X in the time dependent propagator and carrying 
this change all the way through the Green's function treatment. The quantization rule 
and Husimi functions become then 

{S -I){Em) = im+l/2)h (6.66) 

and 

where all quantities are computed with instead of Hi = Ti. 

According to our discussion in section 3, one could also use the Weyl hamiltonian 
and drop I. In this case there would be no e{z) coming from dXjdt. The quantization 
rule in this case becomes 

5(E„) = (m+ l/2)/i , (6.68) 
which is exactly the WKB rule, and the Husimi functions are 

(.) - ^ exp / -(^--^^"^)' 1 (6 69) 

This expression is similar to a formula suggested previously by Kurchan et al. (see 
equations (5.21) and (5.22) in |[Kur89||; see also ||Car92|| for a numerical application). 



The semiclassical energy levels will be slightly different for the three prescriptions, 
eqs. ( |6.62| ), ( |6.66| ), and ( |6.68| ). In fact it can be shown, using arguments similar to those 



in appendix C, that they all coincide up to first order in h. This is in accordance with the 
fact that corrections to the WKB energies are of order h"^ [|Vor77|| . Our quantization rule 



( |6.62| ) does include some corrections of order h"^, but not all of them. These corrections 
might improve the calculation of the energy levels with respect to WKB, especially 
for low lying energy levels. Numerical work is currently being done to explore the 
possibilities. In this we still have one big freedom, the choice of b, the width of the 
coherent state, for each energy level. For the ground state, the variational principle tells 
us that Ti. is always larger than or equal to the true energy. Therefore we must choose 
b so as to minimize the energy. For the other states, however, no such direct rule exists 
and other prescriptions leading to an optimal b have to be devised. 

For the harmonic oscillator, all three formulae give the exact result = (m + 
l/2)hw and the same Husimi distributions, independent of what b is chosen for the 
coherent state width. For the case of 6 = \ffifjj all three semiclassical formulae give 



Pyn{z) = J- exp 

v27r 2 





2 

— m 


- 1/2)^ 


2 


z 


2 



(6.70) 



which is the result we already obtained in subsection 6.2. 
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7. Conclusion 

One inescapable conclusion is that there is more than one semiclassical approximation to 
quantum mechanics in phase space, even when one restricts oneself to the coherent state 
representation. It has often been said that all semiclassical approximations are identical 
in the end, that they always consist in expanding some exponent to second order in the 
vicinity of the classical trajectory and then doing a collection of gaussian integrals, that 
all such expansions should be the same except for the choice of independent variables, 
and therefore that the gaussian integrals should always produce the same result. This 
argument is wrong, the basic reason being that the classical trajectory about which the 
expansion takes place differs from method to method, because the classical hamiltonian 
differs from method to method. The confusion in the literature is in part due to this, 
but not entirely. People have ignored terms that looked small or unfamiliar, even when 
these were clearly part of the approximation they were using. They have done complex 
integrals as if these were real. They have ignored phases. They have been inconsistent 
with the hamiltonian they used, often changing it in midstream. 

Although the number of different, equally valid approximations is actually infinite, 
it is convenient to focus on only three of them, which have been discussed at some length 
in section 3. They correspond to three different choices for the classical hamiltonian 
associated with a given quantal hamiltonian. They are the Weyl hamiltonian Hw, the 
smoothed hamiltonian Hi = H, and the antismoothed hamiltonian H2. It is essential 
to realize that, in all semiclassical approximations, the smoothed hamiltonian is always 
associated with an exponential term containing a special correction to the action which 
we have called I. The temptation to omit this unfamiliar term is great, but it should be 
resisted, as without it one cannot get a correct quantization rule, for instance. Similarly, 
the antismoothed hamiltonian is always associated with —X. The Weyl hamiltonian does 
not come with such a correction term, which in a way makes it the simplest of the three. 
On the other hand, there is no simple approximation involving coherent states which 
yields the Weyl hamiltonian. It occurs in some of the other methods, for instance the 
Wigner-Weyl method mentioned in section 1, and also the Heller method of section 4.3. 
The advantage of the smoothed hamiltonian is that, since it is smoothed, approximations 
based on power series expansions are especially good for it. This is the hamiltonian that 
we have used in most of our work. The antismoothed hamiltonian, lacking both of the 
characteristics that make the other two desirable, does not seem to have been used in 
practice by anyone. 

The first test cases for a semiclassical approximation should be the free particle 
and the harmonic oscillator. All three of the above approximations become exact then. 
Once again, this is so only when the I term is appropriately included. Without I, both 
Hi and H2 are wrong; only Hw remains. Beyond these two systems, the question of 
which of the three approximations gives better results has no single answer; it depends 
very much on the specific quantity calculated. We are planning to address this point in 
another paper. It is true that the three versions differ from each other only by terms of 
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order h, but a difference of tliis order can be essential in the quantization rule, and it 
may also be significant for the long time behavior. 

We shall now summarize our other results. They can be grouped in three 
categories: coherent state propagator, initial value representation (IVR), and energy 
Green's function. We gave in section 2 a very complete derivation of the semiclassical 
coherent state propagator in one dimension, smoothed hamiltonian version. This is eq. 
( p.59| ) or eq. (|2.68|) . Then in section 3 we compared it with several other semiclassical 
approximations. In particular, we obtained an explicit expression for the antismoothed 
hamiltonian version. We believe both of these results to be new. 

In an attempt to provide a rigorous derivation of the well known Herman-Kluk 
formula, we derived in section 4 our own new IVR. Like HK, we set out to calculate 
semiclassically the mixed propagator {x\K{t)\z') in terms of real trajectories starting 
from g(0) = q', p{0) = p'. Our result, eq. ( |4.29| ), is very different from HK's. Also in 
section 4, we considered Heller's old idea for an IVR, which is not based on the coherent 
state formalism, and we derived an explicit expression for it, valid for any smooth 
hamiltonian function of q and p. We found that Heller's IVR and ours are very similar, 
the difference being that Heller's has the Weyl hamiltonian and no X term. We returned 
to the Herman-Kluk IVR in section 5 and pointed out mistakes made in its derivation. 
The nonconservation of the HK norm is a fatal flaw, we think. The HK expression 
does not follow from the semiclassical limit of the coherent state propagator, but it can 
be derived from an ansatz where one assumes a fixed width for the propagated wave- 
packet and one demands that the resulting propagator recover that of Van Vleck when 
transformed to the representation {x\K{t)\x') . We show, however, that the applicability 
of such a formula is very limited and should fail for long times and/or for chaotic systems. 

Finally, in section 6 we derived the corresponding semiclassical approximations 
for the energy Green's function in the coherent state representation. For each of the 
approximations, we used this Green's function to find a quantization rule, which yields 
the energy levels in terms of classical quantities. We also used it to get a semiclassical 
approximation to the Husimi distribution for each level. 

There are at least two ways that one might go farther in this field. One needs to 
carry out numerical comparisons between the various semiclassical approximations and 
determine the conditions under which one or the other might be preferable. We have 
made a small start in this direction, but we leave the results for a future publication. One 
needs also to extend this work to two dimensions, which will allow applications to more 
interesting problems, in particular problems with partially chaotic classical mechanics. 
There too we have made a significant start; it promises to be a rather complicated field 
and we shall say no more about it here. 
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Appendix A. Calculating the Prefactor by the Determinantal Method 

In this appendix we present an alternative way of performing the multiple gaussian 
integral which occurs in the calculation of the semiclassical propagator. In subsection 
2.3, we started by doing the integrals one by one, and then we derived a recursion relation 
between successive integrals. Here we start with the determinant of the quadratic form, 
and we derive recursion relations between determinants of successively higher orders. 
This method is very general and can be applied for any number of dimensions. We use it 
here again for an integral in subsection 3.1 concerning the semiclassical approximation 
of a different kind of coherent state propagator. 

We start with eq. ( 2.23 ), where we had to perform a 2(A^— 1) -dimensional gaussian 
integral of the form 

where M is a complex, positive definite L x L matrix. All the work comes in calculating 
the "prefactor", or the square root of the determinant. We rewrite eq. ( |2.23| ) as 



K,{z",t;z',0) = e^^'*''^ / <{ I I K . (A.2 




The boundary conditions (|2.19|) ensure that j varies only from 1 to — 1. We already 
wrote ( in eq. ( |3.12| ). Here we define the abbreviations 

_ i d'H,{z^^„z,) _ i d'H,{zlz,_,) _ i d'H,{z^_,„z,) 

' h dz^ ' h dzf ' h dz*dz, ^ ^ 



and we write 

N-1 N-1 N-2 N-1 



j=l j=l j=l j=l 



The subscript 1 in eqs. (|A.2| ) and (|A.3|) is the notation of section 3 and means that 



the first form of path integral is being used. We shall calculate K2{z" ,t; z' ,0) as 
well later in this appendix. We use r] (without subscript!) to denote a double- 
sized vector containing all the rj/s as well as their complex conjugates in the order 
Vn-1, V*n-v Vn-2, ■■■,Vi,Vi ■ Then eq.(^) can be rewritten 



C = . ^ (A.5) 
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where 



^(7V-l) .__ 



1 







1 














tBn-i 













tCn-2 - 1 


tAn-2 


1 











1 


rBN-2 


tCn~3 - 1 











tCn~3 - 1 


tAn-s 


1 











1 





(A.6) 

is a tridiagonal band matrix whose dimensionahty is twice the value of its superscript. 

The integration formula ( [A.l[ ) cannot be applied directly to ( |A.2| ), because in the 
latter the variables and the paths are complex. However, this is a problem that we 
already considered and solved at the beginning of subsection 2.3. If we compare carefully 
eq. (|2.25|) with eq. (|2.24| ), we find in the complex case that, for each pair of variables, we 
must introduce an additional minus sign in front of the determinant in the square root. 
It is also clear that the 27r's will cancel out. Hence the correct formula is 



v/(-l)^-MetG(^-i) 



Then the calculation of the determinant proceeds in the following way. 
determinant is expanded with respect to its first column 



(A.7) 
First, the 



detG'(^-i) = rA 



Af-l 



tBn-1 



'tCn-2 







- 1 



tCn-2 — 
'rAN-2 
1 





1 





1 








tAn-3 

1 






1 

tBn- 



1 














'n-2 — 1 


rAN-2 


1 











1 


rBN-2 


tCn-3 - 1 











tCn-3 - 1 


tAn-s 


1 











1 


tBn-3 



(A.8) 

The second determinant can easily be expanded with respect to its first line: it is simply 
det G*^^"^-*. The first determinant, on the other hand, is that of a different quadratic 
form which we call Thus we have the recursion relation 

det G(^-i) = rA^_i det F^^-^) - det G^^-^) _ ^^.9) 

1 T7^Kv.n,„v.,r onno 



Calculating the Prefactor by the Determinantal Method 



70 



Now we can expand detF^^ analogously, which leads to a second recursion relation 

det F(^-i) = tBn^, det G^^-^) - {rCN-2 - 1)' det F^^'^) 

= r5^_i det + 2rC^_2 det F^^-^) _ ^(a.-2) ^ q(^2)_ ^^^^q^ 

Since eventually we shall go to the limit r — > 0, — >■ oo, we may drop the terms of 
order r^. We can get rid of many minus signs by making the definitions 



\N-1 



\N-1 



gN-1 = [-)■■ ^detG(^-^), /^_i = 
Eqs. ( |A.9| ) and (|A.10|) can now be rewritten as 

gN~l = gN-2 + T^N~lfN-l 

In-I = fN-2 — T-BiV-l5'Ar-2 — 2rCAf_2/Af-2 



detF(^-i). 



(A.ll) 



(A.12) 



In the continuous time limit, these recursion relations turn into two first-order 
differential equations 



9 = Af 



f = -Bg- 2Cf . 



(A.13) 



By letting r go to in the determinants themselves, one sees easily that the initial 
conditions are 



giO) = 1 /(O) = . 
One way to solve the differential equations is to introduce the functions 



g = gexp 
Eqs. ( |A.13| ) become now 



C{t')dt' 



f = fexp 



C{t')dt' 



~g = Af + C~g 



f = -B~g-Cf 



(A.14) 



(A.15) 



(A.16) 



with initial conditions g = 1 and / = 0. Replacing the abbreviations A, B, C by 
their definitions in eq. ( A. 3 ), we see that ( |A.16 ) is formally identical to the differential 
equations (|2.47|) if we let g be 5v and f he 5u . The initial values are (5m(0) = 5u' = 
and 5v{0) = 5v' = 1 . 

We need to know g(t), which is the quantity under the square root in ( |A.7|) ; we do 
not need f(t). Since g{t) = 6v", we can use eq. ( p.58| ) with 6v' = 1 to write 

-1 



g(t) = Sv" exp 



C{t')dt' 



i d'S 



exp 



C{t')dt' 



hdu'dv" 

Substituting back into eq. ( |A.7| ) and replacing C by its definition gives 



(A.17) 



Ki{z", t; z' , 0) = e 



h du'dv" 



exp <; - / dt 



dudv j 



(A.18) 
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Finally, by calculating exp{f{z*, z)} in the same way as in section 2 (see eq. ( p. 41 ) and 
subsection 2.5), we obtain a result identical to eq. ( |2.59|) : 



Kiiz",t; z',0) 



1 r..^i 



hdu'dv" ^ ^2h J ^ ^ dudv j 
exp \ ^S{v",u,t) - ^{Wl"^ + \u\^) \ . 



(A.19) 



Let us now repeat this calculation for Klauder and Skagerstam's "second form of 
the path integral" discussed in subsection 3.1. We follow very closely what we just did 
for Ki , with small but important differences. According to eq. (|3.11|) , the quadratic 
form for K2 is 



N-2 



N-l 



N-l 



N-l 



-2C = 2 ^ - 2 J2 vUVj + E ^^^-^1 + E ^^^-^f + 2 E ^'^^'^h (A-20) 



i=i j=i 
with the (different) definitions 

id^H2{z*,Zj] 



i=i 



A, 



h 



dz] 



^ ^ i d'^H2{z*,Zj] 



C, = -r 



id^H2{z*,Zj] 



h dzjdzj 



(A.21; 



The matrix ^ becomes 

/ rAjv-i 1 + tCn-1 

1 + tCn-1 tBn-i 

-1 































V 





-1 

r v4jv_2 1 + tCm~2 

1 + tCn-2 tBn-2 -1 

-1 tAn-3 I + tCn-3 

l+rCjv-3 tBn-s 

(A.22) 



Calling ^ the matrix obtained by removing the first column and the first row of 
G^"^, and expanding both G^~^ and F^~^ with respect to their first lines, we get 



det = tAn., det F^^'^) - (1 + tCn-i)^ det G^^^^) 

det = ri?^_i det G^^-^) _ ^(^-2) _ 



(A.23) 



Once again we get rid of signs with the help of definitions ( |A.11|) and we take the limit 
r ^ 0, — * 00 to get 



g = Af + 2Cg f = -Eg . 



Finally we define 
g = gexp 



Cit')dt' 



f = fexp 



- I C{t')dt' 




(A.24) 



(A.25) 
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Note that, compared to the analogous definitions ( A.15|) , the present ones have a minus 
sign! We obtain again 

~g = Af + Cg } = -Be, - Cf (A.26) 

which are the same equations as (|A.16|) . Therefore the propagator K2 is given again by 
( |A.18| ) except that, this time, there is a minus sign in front of the X term! 
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Appendix B. Proof of eq. (|3.27| ) 

We start by rewriting eq. ( p.24|) as 



A 



fB.ll 



where 0(a;) := f{x) — ih\ogg{x). To calculate A beyond the stationary phase 
approximation (SPA) we expand (f) about its stationary point X to fourth order: 



(x) = 0(X) + j/^^dx" + + j/'^6x^ + ... 



(B.2) 



where 5x = x — X and 0*^"-* = d"0/da;"(A). Next we change the integration variable to 

3*°/'-'"/' 5x , (B.3) 



10(2) I 



2h 



where a is the phase of 0(2). Substituting (|R2D and (Q) into (^) we get 



A 



2h 

\W\' 



gm/4-ia/2gi0(X) 



10(2) I 



3/2 



1 + 



jn0(^) 



6|0(2)|- 



iir/4-3ia/2,,3_i__iA£ifL„iir-2ia,,4 

9|0(2)|3 



dy 



|0( 



2)1 



,i7r/4-ia/2pi<^(X) 



^|0(2)|^ 



5^;i(0('))2 
24|0(2)|3 ' 



-3ia 



10(2) I 



'^(^^ + 2^ (5e-^-(0(^))2 - 3e-2-|0(2)|0(^)) 



(B.4) 



We now proceed to expand the various terms above in powers of h so as to compare this 
result with SPA. The stationary point A of can be written in terms of the stationary 
point Xq of / as 



q(i) 



Also 



7(2)^(0) 



0(A) = /(xo)-^^log^7(a;o) + /i' 



9 



(0) 



(2) _ f(2) 



/^2) + ihf'^^^xi - ih 



gmgio)_giir 



(B.5) 

(B.6) 
(B.7) 
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where g^""^ = d"'g/dx"'{xQ) and f^"^^ = d"//dx"(xo). From the last expression we find 

|0(2)|-^/2 = |/(2)|-i/2(i + o(;^2)) (B_g) 

and 

Substituting these equations into ( |B.4|) we obtain, after some simphfications, the desired 
result, eq. ( p.27| ), with 

, /(2)^(2) - /(3)^(i) ^ 5(/(3))2 - 3/(2)/(^) 
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Appendix C. Cancelation of first order terms in S + 1 

Our goal here is to show that the quantity {S + X) appearing in the phase of the 
propagator, eq. (|2.68| ), can be written as Sw + O^Ji^), where Sw is the action computed 
with the classical Hamiltonian He, which we assume to be the same as the Weyl symbol 
Hw- In order to do that we must relate H to Hw and write the classical trajectories 
governed by 7i in terms of those governed by H^. 

The relation between Hw and Ti. = Hi is given by the second of eqs. ( |3.23|) . 
Expanding the exponential operator we obtain 

Hwiu,v) = n{u,v) - ^^{u,v) + ... . (C.l) 

The second term on the right is what we have called e{u, v) in section 6 (see eq. ( |6.30| )). 
We assume that Hy^ does not depend exphcitly on h, and that TC may be written as a 
power series in h, starting at hP. The same holds for e{u, v), with the series starting at 
h}. Since the next term in the expansion (|C.1|) is of order h"^, eq. (|C.1|) can be rewritten 

as 

Hw{u, v) = n{u, v) - e{u, v) + 0{ff) . (C.2) 

The inverse equation, where Ti. is written in terms of Hw, is 

niu, v) = Hw{u, v) + tw{u, v) + 0{h^) (C.3) 

where ew{u,v) is defined by eq. ( |6.30D with Hw replacing Ti. 

Let u = u{v" ,u' ,t') and v = v{v",u',t') be the solutions of Hamilton's equations 
( p.35| ) with boundary conditions u{v",u',0) = u', v{y" ,u' ,t) = v". Let also uq = 
uo{v", u' , t') and vq = vq{v", u' , t') be the solutions of the same Hamilton's equations but 
with Ti replaced by the classical Hamiltonian Hw and the same boundary conditions, 
uq{v" ,u' = u' , vo{v",u',t) = v" . Since the difference between the functions Tt and 
Hw is of order h we write 

u{v" , u , t') = uo{v" , u , t') + hui{v" , u , t') + 0{h^) 

v{v",u',t') = vo{v",u',t') + hvi{v",u',t') + 0{!f) . (C.4) 

Notice that t is the total propagation time, whereas t' is used for intermediate instants, 
< t' < t. Due to the boundary conditions we have 

ui{v" , u , 0) = 

Vi{v",u',t) = Q . (C.5) 
For the sake of clarity we rewrite formula ( p.52|) for the action: 



S{v" , u , t) : = dt' 



t 

— (-uf — vu) — 7i{u, V, t') 



yKV' + wV). (C.6) 
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We shall now expand the various terms in S{v",u',t) in powers of h. We start with 
TC{u,v) calculated along its own trajectory: 

Hiu, v) = Hw{uo + hui, vo + hvi) + tw{uQ, vq) + 

ou ov 
= Hwiuo, Vo) - ih^VoUi + ih^iioVi + e{uo, Vq) + 0{f\^) (C.7) 

where the dot means d/dt'. Using eqs. ( |U.4| ) we also find 

— {UV - VU) = —{UoVo - VqUo) + —{UlVo + UqVi - VqUi - ViUo) + 0{h^) 

-{iiovo - vqUo) - ih^ivoui - iioVi) + ^^(^oMi - ^^i^o) + • 

(C.8) 



Finally, using eqs. 

ih 



we have 



-{u"v" + u'v') 



|« + K)^" + ju'iv'o + hv[) + Oih') 



(C.9) 



Substituting eqs. (|C.7| )- ([(T9D in eq. ( |C.6|) we see that the second and third terms in eq. 
( p.7| ) cancel the terms in the second parenthesis of ( |C8| ). We get 



S{v",u',t) = Sw{v",u,t) + / dt' 



ih^ d 



-{vqUi - viuo) - e{u,v) 



2 dt 

u'[v" + u'v[) + 0{h^ 



(C.IO) 



Using eq. ( |C.5| ) once again, we cancel the first two terms on the second line against those 
under the total derivative sign. Recalling the definition ( |2.61| ) of X, we finally get 



S{v\ u', t) + I{v\ u', t) = Sw{v", u', t) + 0{h^) . 



(C.ll) 
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